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Preface 


This  user’s  guide  describes  the  computer  program  CBEAMC  which  can  be 
used  for  analysis  of  general  beam-column  structures  supported  and/or  loaded 
by  components  which  interact  with  the  displacements  of  the  beam -column. 

The  woric  in  writing  the  computer  program  and  user’s  guide  was  accomplished 
with  funds  provided  to  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  Vicksburg,  MS,  by  Headquarters,  U.S.  Army  Corps  of  Engineers,  Civil 
Works  Directorate,  imder  the  Structural  Engineering  Research  Program  work 
unit  of  the  Computer-Aided  Structural  Engineering  (CASE)  project. 

The  program  and  user’s  guide  were  written  by  Dr.  William  P.  Dawkins, 
P.E.,  Oklahoma  State  University,  Stillwater,  OK,  under  IPA  agreement  No. 
93-15-M  with  WES. 

The  work  was  managed,  coordinated,  and  monitored  in  the  Information 
Technology  Laboratory  (ITL),  WES,  by  Dr.  Reed  L.  Mosher,  Acting  Chief 
Computer-Aided  Engineering  Division  (CAED).  Mr.  H.  Wayne  Jones,  Chief, 
Scientific  and  Engineering  Applications  Center,  CAED,  is  Project  Manager  for 
the  CASE  Project.  Dr.  N.  Radakrishnan  is  Director,  ITL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 

inches 

2.54 

centimeters 

kips  (force) 

4448.22 

Newtons 

pounds  (force) 

4.4482 

Newtons 

pounds  (force)  per  square  inch 

6894.76 

Pascals 

pounds  (force)  per  square  foot 

47.8803 

Pascals 

1  Introduction 


General 


This  report  documents  a  computer  program-CBEAMC-for  analysis  of 
general  beam-column  structures  supported  and/or  loaded  by  components  which 
interact  with  displacements  of  the  beam  and/or  column. 


Organization  of  Report 

The  remainder  of  this  report  is  divided  into  the  following  chapters: 

a.  Chapter  2.  Describes  the  general  beam-column  system  considered  and 
the  mathematical  model  used  for  analysis. 

b.  Chapter  3.  Presents  the  force-displacement  relationships  for  the  mathe¬ 
matical  model  and  describes  the  computational  procedure  used  for 
solution. 

c.  Chapter  4.  Provides  example  solution  obtained  with  the  program. 


Disclaimer 

The  computer  program  described  in  this  report  has  been  checked  to  ensure 
that  the  results  are  accurate  within  the  limitations  of  the  procedures  employed. 
However,  there  may  be  unusual  situations  which  were  not  anticipated,  and 
these  situations  may  cause  the  program  to  produce  questionable  results.  It  is 
the  responsibility  of  the  user  to  judge  the  validity  of  the  results.  No  responsi¬ 
bility  is  assumed  by  the  author  for  the  design  or  behavior  of  any  structure 
based  on  results  obtained  with  the  program. 
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2  Beam-Column  System 


General 

The  general  beam-column  system  considered  for  analysis  is  shown  in  Fig¬ 
ure  1.  Characteristics  and  assumptions  employed  for  each  of  the  system  com¬ 
ponents  are  described  in  the  following  paragraphs. 


Beam 

The  following  assumptions  of  conventional  beam  bending  theory  are 
employed: 

a.  The  beam  is  composed  of  a  linearly  elastic,  piecewise,  homogeneous 
material  with  modulus  of  elasticity  E. 

b.  The  beam  is  essentially  piecewise  prismatic,  i.e.,  the  cross  section  of 
area  A  and  moment  of  inertia  1  may  vary  along  the  length  but  do  not 
alter  the  basic  assumption  that  plane  cross  sections  remain  plane  after 
loading. 

c.  All  cross  sections  share  an  initially  straight,  common  centroidal  axis; 
the  x-axis  is  shown  in  Figure  1. 

d.  All  cross  sections  gave  a  principal  axis  parallel  to  the  y-axis.  Figure  1. 

e.  Deformations  due  to  shearing  strains  are  negligibly  small. 

/.  Stress  concentrations  at  abrupt  changes  in  cross  sections,  at  concen¬ 
trated  loads  or  springs,  and  at  supports  are  negligibly  small. 
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Supports 


Distributed  Lateral  (y)  Distributed  Axial  (x) 

Spring  Supports  Spring  Supports 


Figure  1 .  General  beam-column  system 


Displacements 

The  following  assumptions  regarding  displacement  of  the  system  are 
employed: 

a.  Displacement  of  points  on  the  beam  centroidal  axis  are  completely 
defined  by  translations  u  (parallel  to  the  x-axis),  and  v  (parallel  to  the 
y-axis),  and  rotation  0  about  and  axis  perpendicular  to  the  x-y  plane. 

b.  All  displacements  are  small  and  do  not  alter  the  basic  geometry  of  the 
system. 

c.  Rotations  0  are  small  enough  to  permit  the  approximations  cos  0=1 
and  sin  0  =  0. 


Applied  Loads 

The  following  assumptions  regarding  external  applied  loads  are  employed: 

a.  Applied  loads  are  unaffected  either  in  magnitude,  direction,  or  point  of 
application  by  the  deflections  of  the  system. 


Chapter  2  Beam-Column  System 


b.  Concentrated  applied  loads  are  assumed  to  be  applied  as  components 
parallel  to  the  x-  and  y-directions  or  as  couples  acting  about  axes  per¬ 
pendicular  to  the  x-y  plane. 

c.  Distributed  applied  loads  are  assumed  to  act  parallel  to  the 
X  and  y  axes. 


Fixed  Supports 

A  fixed  support  is  assumed  to  be  an  external  influence  which  results  in  a 
specific  value  (either  zero  or  nonzero)  of  one  or  more  displacement  compo¬ 
nents  at  a  point  on  the  structure,  regardless  of  other  effects. 


Concentrated  Linear  Spring  Supports 

Concentrated  linear  springs  produce  forces  on  the  stmcture  which  are 
directly  proportional  to  the  displacements  of  the  point  of  attachment.  The  fol¬ 
lowing  assumptions  are  employed: 

a.  Concentrated  linear  translational  springs  may  be  inclined  at  an  angle  (a, 
Figure  1)  to  the  axis  of  the  beam.  The  resisting  force  developed  in  the 
spring  is  directly  proportional  to  the  translation  of  the  point  of  attach¬ 
ment  parallel  to  the  line  of  action  of  the  spring. 

b.  A  concentrated  linear  rotational  spring  applies  a  resisting  couple  pro¬ 
portional  to  the  rotation  of  the  point  of  attachment. 


Distributed  Linear  Spring  Supports 

Distributed  linear  spring  supports  produce  distributed  resisting  forces 
directly  proportional  to  the  displacements  of  the  structure.  The  following 
assumptions  are  employed: 

a.  Distributed  linear  springs  resist  only  lateral  displacements  v  or  axial 
displacements  u. 

b.  The  magnitude  of  the  distributed  resisting  force  depends  only  on  the 
lateral  displacement  v  or  axial  displacement  u  at  any  point,  i.e.,  the 
Winkler  hypothesis. 
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Concentrated  Nonlinear  Spring  Supports 

Concentrated  nonlinear  spring  supports  produce  forces  on  the  structure 
which  depend  on,  but  are  not  directly  proportional  to,  the  displacements  of  the 
point  of  attachment  The  following  assumptions  are  employed: 

a.  Concentrated  nonlinear  translational  springs  may  be  inclined  at  an  angle 
(a.  Figure  1)  with  the  axis  of  the  beam.  The  resisting  force  developed 
in  the  spring  is  a  function  of  the  translation  of  the  point  of  attachment 
parallel  to  the  line  of  action  of  the  spring. 

b.  Nonlinear  rotational  springs  are  not  considered. 

Distributed  Nonlinear  Spring  Supports 

Distributed  nonlinear  spring  supports  produce  distributed  resisting  forces 
which  depend  on,  but  are  not  directly  proportional  to,  the  displacements  of  the 
structure.  The  following  assumptions  are  employed: 

a.  Distributed  nonlinear  springs  resist  only  lateral  displacements  v  or  axial 
displacements  u. 

b.  The  magnitude  of  the  distributed  resisting  force  depends  only  on  the 
lateral  displacement  v  or  axial  displacement  u  at  any  point,  i.e.,  the 
Winkler  hypothesis. 


Characteristics  of  Nonlinear  Springs 

The  force-deformation  relationship  for  a  nonlinear  spring  (either  concen¬ 
trated  or  distributed)  is  assumed  to  be  provided  as  a  piecewise  linear  curve,  as 
shown  in  Figure  2.  For  any  deformation  A,  parallel  to  the  spring  line  of 
action,  the  total  resisting  force  (concentrated  or  distributed)  may  be  expressed 
as 


Q  {or  q)  =  {or  q^)  +  k  ■  A 


(1) 


Hence,  although  an  iterative  solution  is  required,  during  any  iteration  a  non¬ 
linear  spring  may  be  replaced  by  a  combination  of  a  fixed  load  or  q^)  and 
a  linear  spring  with  stiffiiess  k. 
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3  Finite  Element  Model 


General 

Solutions  for  the  general  beam-column  system  are  obtained  with  the  finite 
element  model  described  below. 


Nodes 


Because  the  centroidal  axis  of  the  beam  is  initially  straight,  the  location  of 
a  node  is  completely  defined  by  its  x-coordinate.  For  identification  in  the 
following  paragraphs,  nodes  are  assumed  to  be  numbered  sequentially  starting 
with  node  1  at  the  left  end  of  the  beam,  as  shown  in  Figure  3. 


x-Coord.  ’'3  *i.1  ’'i  ’‘j+1 

Node  No.  1  2  3  .  i-1  i  i+1 

Elem.  No.  1  2  .  i-1  i  . 

Figure  3.  Finite  element  model 


A  minimum  number  of  nodes  are  dictated  by  the  characteristics  of  the 
system: 
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a.  At  the  left  and  right  ends  of  the  beam. 

b.  At  points  of  change  in  beam  properties  (i.e.,  E,  I,  and  /or  A). 

c.  At  the  point  of  application  of  concentrated  loads. 

d.  At  fixed  supports. 

e.  At  the  point  of  application  of  concentrated  spring  supports. 

/.  At  the  beginning  and  end  of  distributed  (either  axial  or  lateral)  loads. 

g.  At  the  beginning  and  end  of  distributed  spring  (either  axial  or  lateral) 
supports. 

Displacements,  internal  forces,  and  other  effects  are  determined  only  at  the 
nodes.  To  provide  more  detailed  description  of  the  variation  of  effects  along 
the  length,  nodes  may  be  defined  at  any  other  location  on  the  beam.  In 
addition,  the  accuracy  of  a  solution  may  be  affected  by  the  number  of  nodes. 
This  aspect  will  be  discussed  subsequently. 


Variations  in  System  Properties 

The  following  assumptions  are  employed  for  system  properties  which  vary 
along  the  length  of  the  structure. 

a.  Properties  E,  A,  and  I  of  beam  cross  section  are  assumed  to  be  constant 
between  adjacent  nodes.  In  effect,  a  tapered  beam  is  replaced  by  a 
piecewise  prismatic  (stepped)  member.  As  discussed  later,  this  replace¬ 
ment  is  accomplished  automatically. 

b.  Distributed  loads  (either  applied  or  resulting  from  the  fixed  load  portion 
of  distributed  nonlinear  springs)  are  assumed  to  vary  linearly  between 
adjacent  nodes. 

c.  Distributed  spring  stiffnesses  (either  distributed  linear  springs  or  the 
linear  stiffness  component  of  distributed  nonlinear  springs)  are 
assumed  to  vary  linearly  between  adjacent  nodes. 


Elements 


An  element  is  defined  as  the  portion  of  the  beam  between  adjacent  nodes. 
For  reference  in  the  following  paragraphs,  an  element  is  identified  by  the  node 
number  at  its  left  end,  as  shown  in  Figure  3. 
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Under  the  assumptions  stated  above,  each  element  is  characterized  by  mod¬ 
ulus  of  elasticity  E,  properties  A  and  I  of  the  cross  section,  and  length  h  and  is 
subjected  to  distributed  loads  and  springs,  as  shown  in  Figure  4. 


Figure  4.  Typical  element 


Interaction  of  internal  axial  stress  resultants  with  lateral  deflections  will 
affect  the  bending  resistance  of  the  element.  When  this  effect  is  included,  the 
bending  resistance  of  the  element  is  developed  using  a  constant  average  axial 
stress  resultant  given  by 

^  "  'X 


The  relationship  between  element  end  forces  and  end  nodal  displacements, 
obtained  by  procedures  given  by  Clough  and  Penzien  (1975),  Przemieniecki 
(1968),  and  Zienkiewicz  (1971),  is 


(3) 


where 
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L  \fx,i  fyji  ^/+l] 

=  (6x1)  vector  of  element  end  forces 

5^  =  (6x)  axial  and  beam  bending  stiffness  matrix  including 
effects  of  axial  stress  resultant  on  bending 

5]^  =  (6x6)  stiffness  matrix  representing  effects  of  distributed  springs 
u'  =  [m-  V,.  e-  v- 


=  (6x1)  vector  end  of  nodal  displacements,  and 

ft  _\f‘  ft  i  ft  ft  i 

J-e  Jexj  ^ey,i  '^ej  ^ex,i*l  ■'ey,i*l  ^e,/+l 

=  (6x1)  vector  of  fixed  end  forces  due  to  distributed  loads. 


Coefficients  for  matrices  and  ^are  given  in  Figures  5,  6,  and  7 

and  Table  1. 


Node  Equilibrium 


A  typical  node,  shown  in  Figure  8,  is  subjected  to  : 

a.  Concentrated  forces  Fy and  C,,  where  F^  ■  and  Fyj  include  both 
applied  external  concentrated  loads  and  the  fixed  force  component  of 
all  concentrated  nonlinear  springs  attached  to  the  node. 

b.  One  or  more  concentrated  linear  translation  springs  which  include  the 
linear  spring  components  of  concentrated  nonlinear  springs  attached  to 
the  node. 

c.  One  or  more  linear  rotation  springs. 

d.  Element  end  forces  finm  elements  on  either  side  of  the  node. 


Equilibrium  of  node  i  is  expressed  by 
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S- 

h 


a.  Explicit 


S'  = 


S' 

S' 

2e.i+I.i 


a^- 


SYM 


]_ 

h 


asl 


-A 

0 

0 

A 


S' 

^E.i.i+l 

S' 

^E.i+l,i+l 


b.  Partitioned  symbolic 


0 


I 

a4l 

0 

I 

asl 


where:  E  =  modulus  of  elasticity 

A  =  area  of  cross  section 
I  =  moment  of  inertia  of 
cross  section 

a,  -  =  constants  given  in 

Table  1. 


Figure  5.  Element  axial  and  bending  beam  stiffness  matrix 


Figure  6.  Contribution  of  distributed  springs  to  element  stiffness  matrix 
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-(20q,.+10q,,,,,)' 

(bi^y.i +t>2qy.i-^|)  Where  b,,  b,,  by,  b4 

^(^sQy.i  +  b4Qy.i+i )  are  constants  given  in 
“(lOq^.j  +20q^  ,^,)  ^ 

(b2qy,+b,qy,„) 

-h(b4q^,,+b3qy,„) 


a.  Explicit 


n-  ■ 

If 

i-e.i  +  l 


b.  Partitioned  symbolic 


Figure  7.  Element  fixed  end  forces  due  to  distributed  loads 


k^jcos^a  A:^,/ sin  a  cos  a  q 

/y./^  j  /j),/ 1  ■*■52  ^4,,/ sin  a  cos  a  k^  -sii?a  0 


M  F 

/  X,l 


h  r>..' 


e,  C; 


where  the  summation  X  extends  over  all  concentrated  springs  attached  to 
node  i. 
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Table  1 

Constants  for  Element  Stiffness  Matrix  and  Fixed  End  Forces 
(P  =  Axial  Internal  Stress  Resultant) 


Coefficient 

P  =  0 

P  =  Compression 

P  =  Tension 

a, 

12 

<b's/A 

d)'s/A 

a2 

6 

<b"(l-c)/A 

-d)"(l-c)/A 

aj 

4 

<b^(s-Oc)/ A 

-^^(s-Oc)/ A 

a^ 

2 

0^(0 -s)/ A 

-<b'('b-s)/A 

b, 

-21 

10[-3(4  +  O')(l-c) 
+Os(6-t-2<b^)]/A 

10[3(4-(I)')(l-c) 
+08(6-2(1)^)] /A 

b2 

-9 

10[3(4-d)-)(l-c) 

-cbs(6-0-)]/A 

10[-3(4  +  (I)^)(l-c) 

-08(6  +  0')] /A 

bs 

-3 

10[9<I)s-O-(l  +  2c) 
-12(l-c)]/A 

10[9O8-O'(1  +  2c) 
+12(l-c)]/A 

b4 

-2 

10[30s-0-(2  +  c)]/A 

10[30s-0'(2  +  c)]/A 

A 

... 

0^(2-2c-<bs) 

O' (2 -2c +  08) 

(b 

— 

VPhVEI 

VPh'/EI 

s 

... 

sinO 

8inhO 

c 

... 

cosd> 

coshO 

Substitution  for  end  force  vectors and//  in  terms  of  element  stiffness 
matrices  using  Equation  3  and  Figures  5,  6,  and  7  results  in  the  following 
governing  equation  which  must  be  satisfied  at  every  node. 
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Concentrated  Forces 


Figure  8.  Typical  node 


Node  Displacements 


*  til/  "  ilj  *  4/./  *  4./  * 

*i4i,/-i  *  4.1 1  ■‘4.1  -4,  41  4./ 

Evaluation  of  Equation  5  at  every  node  in  the  model  results  in  a  system  of 
simultaneous  equations  (3n  equations  for  a  model  with  n  nodes)  which  are 
solved  for  the  nodal  displacements.  When  nodal  displacements  are  known, 
internal  forces,  spring  forces,  and  reactions  may  be  determined  for  the  system. 

Iteration 


A  single  solution  of  the  simultaneous  equations  produces  final  results  for 
systems  possessing  the  following  characteristics: 

a.  Only  linear  spring  supports  are  present;  and 
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b.  The  effect  of  axial  stress  resultant  on  element  bending  stiffness  is 
excluded:  or 

c.  The  effect  of  axial  stress  resultant  on  element  bending  stifftiess  is 
included,  axial  distributed  springs  are  absent,  and  all  concentrated  trans¬ 
lational  springs  are  perpendicular  to  the  beam  axis. 

All  other  systems  require  an  iterative  solution  to  evaluate  axial  stress  resul¬ 
tants  and/or  fixed  force-linear  spring  components  for  nonlinear  springs.  The 
iterative  solution  process  is  initiated  by  evaluating  axial  stress  resultants  and 
nonlinear  spring  characteristics  for  zero  displacements.  The  simultaneous 
equations  are  evaluated  and  solved  for  a  new  estimate  of  system  displacements. 
On  each  iteration,  the  solution  of  the  preceding  iteration  is  used  for  evaluation 
of  system  properties.  This  process  is  continued  until  the  results  of  two  succes¬ 
sive  iterations  differ  by  an  acceptably  small  amount. 


Effect  of  Node  Spacing  on  Solution 

The  finite  element  procedure  described  produces  “exact”  solutions  for  sys¬ 
tems  possessing  the  following  characteristics: 

a.  The  beam  is  piecewise  prismatic. 

b.  Only  concentrated  external  loads  and  linearly  varying  lateral  (y)  loads 
are  present.  Or,  if  the  effect  of  axial  stress  resultant  on  element  bend¬ 
ing  stiffness  is  excluded,  distributed  axial  (x)  loads  vary  linearly. 

c.  Only  concentrated  springs  are  present  and  ,  if  the  effect  of  axial  stress 
resultant  on  element  bending  stiffness  is  included,  all  concentrated 
translation  springs  are  peipendicular  to  the  beam  axis.  (Note:  The 
accuracy  of  a  solution  including  the  effect  of  axial  stress  resultant  in 
the  presence  of  inclined  translational  springs  may  be  affected  by  the 
iteration  convergence  tolerance  but  not  by  node  spacing.) 

d.  Adjacent  elements  are  approximately  equal  in  length.  Adjacent  ele¬ 
ments  having  drastically  different  lengths  may  result  in  significant 
round-off  errors  in  the  solution  of  simultaneous  equations. 

In  all  other  systems,  the  number  of  nodes  (and  elements)  used  in  the  finite 
element  model  affects  the  accuracy  of  the  solution.  In  general,  as  the  number 
of  nodes  (and  elements)  is  increased,  the  solution  tends  to  converge  to  the 
“exact”  solution.  There  is  no  “rule-of-thumb”  to  provide  the  necessary  number 
of  nodes  for  acceptable  results.  It  may  be  necessary  to  perform  several  solu¬ 
tions  for  various  numbers  of  nodes  and  node  spacings  to  ensure  that  an 
adequate  solution  has  been  obtained. 
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4  Computer  Program 


General 

A  computer  program  implementing  the  analytical  process  described  in  the 
previous  chapter  has  been  written  in  the  FORTRAN  programming  language 
with  all  arithmetic  operations  performed  in  double  precision. 


Input  Data 

The  program  has  been  written  to  generate  automatically  intermediate  data 
values  from  a  minimum  of  user  input  data.  Details  of  required  user  input  data 
are  presented  in  the  Input  Guide,  Appendix  A.  The  processes  and  assumptions 
used  in  converting  user  input  data  to  intermediate  data  values  are  described 
below. 


Global  Coordinate  System 

A  horizontal  beam  centroidal  axis  (the  x-axis,  positive  to  the  right)  is 
assumed.  The  global  y-axis  is  positive  upward.  The  origin  of  the  x-y  coordi¬ 
nate  system  is  arbitrary.  A  local  coordinate  system  associated  with  nonlinear 
springs  is  described  later. 


Displacement  and  Load  Sign  Conventions 

Positive  directions  for  displacements  and  loads  are  as  follows: 

a.  Nodal  displacements  u  and  v  are  positive  if  the  node  translates  in  the 
positive  X-  and  y-directions,  respectively. 

b.  Nodal  rotation  is  positive  if  the  node  rotates  counterclockwise. 
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c.  External  loads  (concentrated  Fj^  and  Fy;  or  distributed  and  Qy)  are 
positive  if  they  act  in  the  positive  x-  and  y-directions,  respectively. 

d.  A  concentrated  applied  couple  (C)  is  positive  if  it  tends  to  rotate  the 
node  counterclockwise. 

e.  Concentrated  or  distributed  linear  spring  stiffnesses  are  always  positive. 


Nonlinear  Spring  Conventions 

As  stated  previously,  the  characteristics  of  nonlinear  springs  are  assumed  to 
be  described  by  a  piecewise  linear  curve  giving  spring  resisting  force  as  a 
function  of  spring  deformation.  The  sign  conventions  assumed  for  each  of  the 
three  types  of  nonlinear  springs  permitted  by  the  program  are  as  follows. 


Concentrated  Nonlinear  Springs 

Concentrated  nonlinear  springs  are  assumed  to  resist  nodal  translation  dis¬ 
placement  components  parallel  to  the  spring  line  of  action.  Positive  resisting 
forces  and  spring  deformations  for  concentrated  springs  are  defined  for  a  local 
x’-coordinate,  as  shown  in  Figure  9.  The  origin  of  the  local  x’-coordinate  is  at 
the  point  of  attachment  (node);  the  positive  x’-direction  extends  along  the 
spring  line  of  action  away  from  the  point  of  attachment  The  inclination  of  the 
spring  is  given  by  the  angle  a  between  the  global  x-axis  and  the  local  x’-axis 
(positive  counterclockwise). 

Spring  resisting  forces  are  positive  if  they  act  in  the  positive  local  x’-direc- 
tion.  To  illustrate,  consider  a  spring  in  an  initial  state  of  compression  at  zero 
nodal  displacement  The  spring  therefore  would  produce  a  force  acting  on  the 
node  in  the  negative  x’-direction  (i.e.,  negative  force),  -Fj,  at  zero  spring 
deformation,  Aj  =  0,  as  illustrated  by  point  (1)  in  Figure  9b.  If  the  point  of 
attachment  undergoes  nodal  displacements  u  and  v  to  produce  a  positive  defor¬ 
mation  +A2,  Figure  9a,  the  compression  (i.e.,  negative  spring  force),  -F2, 
would  increase  as  illustrated  by  point  (2)  in  Figure  9b.  Further  positive  defor¬ 
mation  to  A3  would  result  in  a  negative  resisting  force,  -F3,  as  shown  by 
point  (3).  If  the  nodal  displacements  produce  deformation  in  the  negative  x’- 
direction,  e.g.,  -A4,  the  compressive  (negative)  force  in  the  spring  would 
reduce  to,  -F4,  point  (4).  At  some  negative  deformation,  -A5,  the  initial  com¬ 
pression  might  be  totally  relieved,  i.e.,  F5  =  0,  as  shown  by  point  (5).  Further 
negative  deformation,  -Ag,  might  produce  tension  in  the  spring  (+Fg),  as  illus¬ 
trated  by  point  (6).  The  coordinate  values.  A,  F,  of  each  point  on  the  deforma¬ 
tion-resistance  curve  are  provided  as  input.  A  linear  variation  between  input 
deformation-resistance  coordinates  is  assumed. 
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a.  Concentrated  nonlinear  spring 


Spring  Resistance 
in  x‘-Direction 


Spring  Deformation 
in  x'-Direction 


b.  Resisting  force-deformation  curve 


Figure  9.  Characteristics  of  concentrated  nonlinear  spring 


Distributed  Nonlinear  Springs 

Distributed  nonlinear  springs  produce  distributed  forces  proportional  to  the 
translation  displacements  of  the  nodes.  The  deformation  of  the  spring  is  equal 
to  and  has  the  same  sign  as  the  nodal  displacement.  The  distributed  force 
produced  on  the  beam  by  a  distributed  nonlinear  spring  is  positive  if  the  force 
acts  in  the  positive  global  x-  or  y-coordinate  directions.  Distributed  spring 
characteristics  are  assumed  to  be  described  by  pairs  resisting  force-deformation 
coordinate  values  similar  to  those  for  concentrated  springs. 
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Data  Generation 


Concentrated  applied  loads,  concentrated  springs,  and  fixed  support  occur  at 
isolated  points  on  the  structure.  All  other  data  (beam  properties,  E,  A,  and  I, 
and  distributed  loads  and  springs)  are  distributed  over  a  range  of  x-coordinates. 
The  user  is  required  to  provide  the  locations  (x-coordinates)  at  which  concen¬ 
trated  effects  occur  and  the  terminal  values  (beginning  and  end)  of  distributed 
quantities.  These  define  the  “dictated”  nodes  in  the  model.  The  program  also 
provides  for  defining  nodes  intermediate  to  the  dictated  nodes.  The  following 
paragraphs  describe  the  procedures  and  assumptions  used  to  obtain  system 
characteristics  at  the  intermediate  nodes. 


Properties  of  Beam  Cross  Section 

Properties  A  and  1  of  beam  cross  section  are  provided  at  the  terminals  of  a 
distribution.  Figure  10  illustrates  for  cross  section  area  the  steps  used  to  con¬ 
vert  input  data  to  element  properties.  The  user  has  specified  values  of  cross 
section  area  Aj  and  Aj^  at  Xj  and  Xj^,  respectively,  resulting  in  “dictated” 
nodes  i  and  i44.  Through  node  spacing  data  (see  below),  intermediate  nodes 
have  been  defined  at  i+1,  i+2,  and  i+3.  The  initial  assumption  is  made  that 
cross  section  area  varies  linearly  from  Aj  at  Xj  to  Aj^,  as  shown  in  Fig¬ 
ure  10b,  to  result  in  values  of  cross  section  areas  Aj^j,  Aj^j*  ^i+3  ^+i’ 

Xj^.2,  and  Xj^3,  respectively.  The  distribution  is  further  converted  to  the  piece- 
wise  prismatic  representation  shown  in  Figure  10c  by  assigning  to  each  ele¬ 
ment  covered  by  the  distribution  a  constmt  area  equal  to  the  average  of  the 
areas  at  each  end  of  the  element,  e.g.,  A*'*"^  =  (Aj^j  +  Note  that  E  is 

required  to  be  constant  over  the  distribution. 

Properties  data  of  beam  cross  section  are  assumed  to  define  the  extreme 
(left  and  right)  ends  of  the  beam-column  system  to  be  analyzed.  Nonzero  val¬ 
ues  of  section  properties  E,  A,  and  1  must  be  supplied  for  every  point  on  the 
structure  between  these  extremes.  Any  other  distributed  or  concentrated  data 
which  fall  beyond  these  extremes  are  ignored  by  the  program. 


Node  Spacing  Data 

Node  spacing  data  provide  for  defining  nodes  intermediate  to  the  “dictated” 
nodes.  The  user  specifies  a  maximum  node  spacing,  to  be  used  in  a 
range  of  x-coordinates.  For  example,  consider  the  range  of  x-coordinates 
between  “dictated”  nodes  Xj  and  x^^  in  Figure  10a.  The  range  of  coordinates 
Xj  to  Xj^  will  be  subdivided  into  elements  of  equal  length  as  follows.  The 
number  of  elements  is  equal  to  the  largest  integer  given  by 
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Dictated  Node 
Start  Distribution 


Dictated  Node 
End  Distribution 


Intermediate  Nodes 


Node  i 


Element 


X-Coord.  X. 


Input  A. 


a.  Notation 


Node  i 


b.  Assumed  linear  distribution 


Element 


c.  Piecewise  prismatic  representation 


Figure  10.  Distribution  of  cross-sectional  properties  of  beam 


n  =  -  Xi)/H^  +  0.9 
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and  the  length  of  elements  in  this  region  will  be 

h  =  (x/+4  -  jc,)//i  O) 

If  is  greater  than  or  equal  to  (x-^^  -  Jt,),  then  the  entire  region  will  be 
trea^  as  a  single  element. 

If  a  “dictated”  node  falls  between  the  limits  of  a  range  of  node  spacing 
data,  intervals  to  the  left  and  right  of  the  dictated  node  are  treated  as  separate 
ranges  for  intermediate  node  spacing. 

Distributed  Loads 

Distributed  loads  vary  linearly  between  the  values  at  the  beginning  and  end 
of  the  distribution  provided  as  input 


Distributed  Linear  Springs 

The  program  provides  for  variations  in  linear  spring  stiffness  coefficients 
over  an  interval  to  X2  given  by  the  general  expression 

k^a  +  b2^  (8) 

where 

k  =  stiffness  coefficient  (force/length/displacement)  of  distributed  x  or  y 
spring  at  any  point  Xj  to  Xj 

a  =  stiffness  at  itj,  always  positive 

b  =  constant  with  units  such  that  bf  has  units  of  stiffness;  b  may  be  posi¬ 
tive  or  negative;  however,  only  positive  values  of  k  are  permitted 

z  =  distance  measured  from  Xj 

c  =  dimensionless  constant 

For  the  input  values  x,-,  x,-^3,  a,  b,  and  c.  the  program  calculates  at  each 
intermediate  node  the  distributed  spring  stiffness  from  the  general  expression 
given  in  Equation  8,  as  illustrated  in  Figure  11. 


Chapter  4  Computer  Program 


Figure  1 1 .  Distributed  linear  spring  stiffnesses 


Distributed  Nonlinear  Springs 

Distributed  nonlinear  springs  are  described  by  resisting  force-deformation 
curves  for  each  end  of  the  distribution,  as  illustrated  in  Figure  12.  As  stated 
previously,  a  nonlinear  spring  is  represented  by  a  combination  of  a  fixed  force 
and  a  linear  spring.  The  components  of  the  nonlinear  spring  at  intermediate 
nodes  are  obtained  as  follows.  Consider  node  /+1  in  Figure  12.  The  spring 
deformation  at  node  r+1,  Aj^j  (either  Uj^.j  or  Vj^j)  is  used  with  the  input  curve 
at  node  i  to  obtain  |,  and  and  with  the  curve  at  node  /+3  to  obtain 
and  The  values  at  node  1+1  are  then  calculated  from  linear  interpolation  as 

•  (^,+3  -  ^1+1 )  <io.r  '  (^,+1  - 

and 

^,+1  +  [^1  '  (^,+3  ■  ^/+l)  ■  (^1+1  ■  ^,)]/^+3  ■ 

Output  Data 

Output  data  are  provided  in  three  parts.  Output  may  be  directed  to  the  user 
terminal,  to  a  data  file,  to  both  or  some  parts  of  the  output  may  be  omitted 
entirely. 
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Echoprint  of  Input  Data 


This  section  contains  a  tabular  listing  of  all  input  data  currently  available  to 
the  program.  This  section  is  output  at  the  option  of  the  user. 


Summary  of  Results 

This  section  contains;  maximum  positive  and  maximum  negative  values  of 
displacements  and  internal  forces  and  the  x-coordinates  at  which  the  maxima 
occur,  reactions  at  fixed  supports;  and  forces  in  concentrated  linear  and/or 
nonlinear  springs.  This  section  is  always  output  to  the  user  terminal,  to  an 
output  file,  or  lx)th. 


Complete  Results 

This  section  contains  a  tabulation  of  displacements,  internal  stress  resul¬ 
tants,  and  forces  in  distributed  springs  at  nodes  contained  within  range  of 
x-coordinates  specified  by  the  user.  This  section  may  be  omitted  entirely. 
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Graphics  Output 


Graphical  displays  of  input  data  consist  of: 

a.  Schematic  of  variation  of  axial  and  bending  stiffnesses  (EA  and  El, 
respectively). 

b.  Schematic  showing  locations  of  fixed  and/or  concentrated  spring 
supports. 

c.  Schematic  showing  locations  and  magnitudes  of  applied  loads. 

d.  Schematic  showing  ranges  of  linear  and/or  nonlinear  distributed  springs. 

Graphical  displays  of  results  include: 

a.  Axial  and  lateral  displacement  diagrams. 

b.  Shear  and  bending  moment  diagrams. 

c.  Forces  in  distributed  springs. 

Sign  Conventions  for  Output 

Sign  conventions  used  for  output  data  are  presented  in  Table  2. 
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Table  2 

Output  Sign  Conventions 

Quantity 

Sign  Convention 

Translations  u  and  v 

Positive  if  node  moves  in  positive  x-  or 
y-coordinate  direction,  respectively 

Rotation 

Positive  if  node  rotates  counterclockwise 

Axial  Internal  Stress  Resultant 

Positive  if  tension 

Shear  Stress  Resultant 

Positive  if  shear  force  tends  to  move  left  end 
of  segment  of  beam  in  positive  y-direction 
(i.e.,  up  on  left,  down  on  right  end  of 
segment) 

Bending  Moment 

Positive  if  produces  compression  in  plus 
y-face  (top)  of  beam 

X-,  Y-Reaction  at  Fixed  Support 

Positive  if  acts  in  positive  x-  or  y-coordinate 
direction,  respectively 

Moment  Reaction  at  Fixed  Support 

Positive  if  acts  counterclockwise  on  beam 

Force  in  Concentrated  Linear  Translation 

Spring 

Positive  if  tension 

Moment  in  Linear  Rotation  Spring 

Positive  if  acts  counterclockwise  on  beam 

Force  in  Concentrated  Nonlinear  Translation 
Spring 

Sign  determined  according  to  input  resisting 
force-deformation  curve  coordinates 

Force  in  Distributed  Linear  Spring 

Positive  if  tension 

Force  in  Distributed  Nonlinear  Spring 

Sign  determined  according  to  input  resisting 
force-deformation  curve  coordinates 
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5  Example  Solutions 


The  examples  presented  below  are  intended  only  to  illustrate  the  use  of  the 
program  and  are  not  to  be  interpreted  as  a  guide  for  application  of  the 
program. 


Example  1 :  Fixed  End  Beam 

The  fixed  end  beam  shown  in  Figure  13  was  analyzed  for  a  succession  of 
loading  conditions  described  in  the  input  data  file  shown  in  Figure  14. 


y 


Cross  Section:  E  =  30- 10*  psi 
A  =  4  sq  in. 

1  =  40  in." 


Figure  13.  Fixed  end  beam  for  Example  1 

a.  The  basic  problem  is  described  in  lines  1000  through  1120  of  the  file. 
During  execution,  effects  of  axial  load  on  bending  were  excluded  for 
the  solution  of  this  part. 
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1000  'EXAMPLE  1  —  FIXED  END  BEAM  ANALYSIS 
1010  'EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 
1020  BEAM 

1030  0  120  30. E6  4  40 
1040  NODES  F 
1050  -10  20  5 
1060  LOADS  F  K 
1070  D  y  0  1.2  10 
1080  C  0  40  0  0 
1090  FIXED  F 
1100  0  FREE  0  0 
1110  10  0  0  0 
1120  FINISH  RERUN 

1130  'EXAMPLE  lA  —  SAME  AS  EXAMPLE  1 

1140  'EXCEPT  INCLUDE  EFFECTS  OF  AXIAL  LOAD  ON  BENDING 
1150  'AXIAL  LOAD  =40  KIPS 
1160  FINISH  RERUN 

1170  'EXAMPLE  IB  —  INCREASE  AXIAL  LOAD  TO  2000  KIPS 
1180  LOADS  ADD 
1190  C  0  1960  0  0 
1200  FINISH  RERUN 

1210  'EXAMPLE  1C  —  INCREASE  AXIAL  LOAD  TO  3000  KIPS 
1220  LOADS  ADD 
1230  C  0  1000  0  0 
1240  FINISH  RERUN 

1250  'EXAMPLE  ID  —  INCREASE  AXIAL  LOAD  TO  3289  KIPS 
1260  'THIS  AXIAL  LOAD  IS  SLIGHTLY  LESS  THAN  BUCKLING  LOAD 
1270  LOADS  ADD 
1280  C  0  289  0  0 
1290  FINISH  RERUN 

1300  'EXAMPLE  IE  —  INCREASE  AXIAL  LOAD  TO  3290  KIPS 

1310  'THIS  AXIAL  LOAD  IS  SLIGHTLY  GREATER  THAN  BUCKLING  LOAD 

1320  LOADS  ADD 

1330  C  0  1  0  0 

1340  FINISH 


Figure  14.  Input  file  for  Example  1 

b.  On  the  first  rerun.  Example  lA,  only  the  problem  heading  was  changed 
to  that  shown  on  lines  1140  through  1160  of  the  data  file.  All  other 
data  remain  unchanged.  During  execution,  the  effects  on  bending  due 
to  the  40-kip^  axial  load  were  included  in  the  solution. 

c.  On  subsequent  reruns,  lines  1170  through  1340,  the  axial  load  was 
increased  by  adding  concentrated  load  data  to  each  preceding  data  set. 


Interactive  control  of  the  program  is  shown  in  Figure  15. 


^  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  vi. 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.35.51 


ARE  INPUT  DATA  TO  BE  READ  FROM  YOUR  TERMINAL  OR  A  FILE? 
ENTER  'TERMINAL'  OR  'FILE', 
f 

ENTER  INPUT  FILE  NAME  (64  CHARACTERS  MAXIMUM) . 
xOOSOdl 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR  TERMINAL, 

TO  A  FILE,  TO  BOTH,  OR  NEITHER? 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER', 
f 

ENTER  OUTPUT  FILE  NAME  (64  CHARACTERS  MAXIMUM). 
x0050dl . out 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  INPUT  DATA? 

ENTER  'YES'  OR  'NO'. 

y 

PLOT  OPTIONS: 


OPTION  PLOT  GENERATED 

1  . AXIAL  AND  BENDING  STIFFNESSES 

2  . FIXED  AND/OR  CONCENTRATED  SPRING  SUPPORTS 

3  . LOADS 

4  . DISTRIBUTED  LINEAR  SPRINGS 

5  . SCHEMATIC  OF  DISTRIBUTED  NONLINEAR  SPRINGS 


ENTER  PLOT  OPTION  (1  TO  5)  OR  'FINISHED'. 

1  (Note:  Plots  generated  by  options  1,  2  and  3  are  shown  in 

Figures  17  through  19.  Items  covered  by  options  4 
and  5  are  not  present  in  this  example.) 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO', 
n 

SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 
X0050D1.0UT 
OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH'. 


Figure  15.  Interactive  program  control  for  Example  1  (Sheet  1  of  7) 

Echoprints  of  input  data  were  requested  only  for  the  first  and  last  in  the 
succession  of  problems.  Figures  16  and  23,  respectively.  Plots  of  input  data 
for  the  initial  run  of  the  sequence  are  shown  in  Figures  17  through  19. 
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DESIRED  OUTPUT  UNITS?  ENTER: 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE', 
d  (Note:  Summary  of  results  is  shown  in  Figure  20.) 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO'.. 

y 

PLOT  OPTIONS: 


OPTION  PLOT  GENERATED 

1  . AXIAL  AND  LATERAL  DISPLACEMENTS 

2  . SHEAR  FORCES  AND  BENDING  MOMENTS 

3  . FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 

4  . FORCES  IN  DISTRIBUTED  NONLINEAR  SPRINGS 


ENTER  PLOT  OPTION  (1  TO  4)  OR  'FINISHED'. 

1  (Note:  Plots  for  options  1  and  2  are  shown  in  Figures 

21  and  22.  Items  covered  by  options  3  and  4 
are  not  present  in  this  example.) 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

ENTER  'YES'  OR  'NO', 
n 

OUTPUT  COMPLETE. 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR  TERMINAL, 

TO  FILE 

X0050D1.0UT 

TO  BOTH  OR  NEITHER? 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER', 
n 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 


Figure  15.  (Sheet  2  of  7) 

The  summary  of  results  was  output  for  each  problem,  Figures  20  and  23. 
Because  the  model  used  to  represent  the  beam  contains  only  three  nodes  and 
two  elements,  no  additional  information  would  have  been  provided  by  a  com¬ 
plete  tabulation. 

Plots  of  the  results  are  shown  in  Figures  21  and  22.  The  plots  of  lateral 
deflection.  Figure  21,  and  bending  moment.  Figure  22,  emphasize  the  fact  that 


Chapter  5  Example  Solutions 


SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 
X0050D1.0UT 
OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH', 
f 

DESIRED  OUTPUT  UNITS?  ENTER: 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE', 
d  (Note:  Summary  of  results  is  shown  in  Figure  20.) 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

ENTER  'YES'  OR  'NO', 
n 

OUTPUT  COMPLETE. 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR  TERMINAL, 
TO  FILE 
X0050D1.0UT 
TO  BOTH  OR  NEITHER? 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER', 
n 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 
X0050D1.0UT 
OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH'. 


Figure  15.  (Sheet  3  of  7) 

results  are  calculated  only  at  the  nodes  of  the  model.  The  use  of  a  larger  num¬ 
ber  of  nodes  would  have  resulted  in  more  appropriately  curved  diagrams. 

The  deflections  and  forces  presented  are  “exact”  solutions  for  each  loading. 
The  effect  of  axial  load  on  bending  is  to  increase  the  deflections  and  bending 
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DESIRED  OUTPUT  UNITS?  ENTER: 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE', 
d  (Note:  Summary  of  results  is  shown  in  Figure  20.) 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

ENTER  'YES'  OR  'NO', 
n 

OUTPUT  COMPLETE. 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHO PRINTED  TO  YOUR  TERMINAL, 
TO  FILE 
X0050D1.0UT 
TO  BOTH  OR  NEITHER? 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER', 
n 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 
X0050D1.0UT 
OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH', 
f 

DESIRED  OUTPUT  UNITS?  ENTER: 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE', 
d  (Note:  Summary  of  results  is  shown  in  Figure  20.) 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO'. 


Figure  15.  (Sheet  4  of  7) 

moments.  As  the  axial  load  approaches  the  theoretical  buckling  load 
(3,289.87  kips),  deflections  and  moments  attain  extremely  large  values. 

However,  for  an  axial  load  below  buckling,  deflections  and  moments  are 
consistent  with  the  direction  of  the  applied  lateral  load.  For  an  axial  load 
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n 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO'. 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 
ENTER  'YES'  OR  'NO'. 

n 

OUTPUT  COMPLETE. 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TO  FILE 

X0050D1.0UT 

TO  BOTH  OR  NEITHER? 

TERMINAL, 

n 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER'. 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO'. 

n 

DO  YOU  WANT  TO  PLOT  INPUT  DATA? 

ENTER  'YES'  OR  'NO'. 

n 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 
ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

STIFFNESS 

y 

SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 

X0050D1.0UT 

OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH'. 

DESIRED  OUTPUT  UNITS?  ENTER: 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE'. 

d 

(Note:  Summary  of  results  is  shown  in 

Figure  20. ) 

n 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO'. 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO'. 

n 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

n 

ENTER  'YES'  OR  'NO'. 

Figure  15.  (Sheet  5  of  7) 


slightly  greater  than  the  buckling  load,  deflections  and  moments  undergo  a 
reversal  in  sign  (i.e.,  are  inconsistent  with  lateral  load  direction)  indicating  that 
buckling  has  occurred. 
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OUTPUT  COMPLETE. 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR  TERMINAL, 

TO  FILE 

X0050D1.0UT 

TO  BOTH  OR  NEITHER? 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER'. 

y 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 
X0050D1.0UT 
OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH', 
f 

DESIRED  OUTPUT  UNITS?  ENTER; 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE', 
d  (Note:  Summary  of  results  is  shown  in  Figure  20.) 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

ENTER  'YES'  OR  'NO', 
n 

OUTPUT  COMPLETE. 

INPUT  COMPLETE. 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR  TERMINAL, 

TO  FILE 

X0050D1.0UT 

TO  BOTH  OR  NEITHER? 

ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER', 
f  (Note:  Echoprint  for  this  part  is  shown  in  Figure  23.) 


Figure  15.  (Sheet  6  of  7) 
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DO  YOU  WANT  TO  PLOT  INPUT  DATA? 


ENTER  'YES'  OR  'NO', 
n 

INPUT  COMPLETE. 

DO  YOU  WANT  TO  CONTINUE  WITH  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS 
INCLUDED  IN  THE  SOLUTION? 

ENTER  'YES'  OR  'NO'. 

y 

SOLUTION  COMPLETE. 

DO  YOU  WANT  RESULTS  PRINTED  TO  YOUR  TERMINAL, 

TO  FILE 
X0050D1.0UT 
OR  BOTH? 

ENTER  'TERMINAL',  'FILE',  OR  'BOTH', 
f 

DESIRED  OUTPUT  UNITS?  ENTER: 

'DEFAULT',  'CONSISTENT',  OR  'SELECTIVE', 
d  (Note:  Summary  of  results  is  shown  in  Figure  23.) 

DO  YOU  WANT  COMPLETE  RESULTS  OUTPUT? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  TO  PLOT  RESULTS? 

ENTER  'YES'  OR  'NO', 
n 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

ENTER  'YES'  OR  'NO', 
n 

OUTPUT  COMPLETE. 

DO  YOU  WANT  TO  EDIT  INPUT  DATA? 

ENTER  'YES'  OR  'NO', 
n 

LAST  INPUT  FILE  PROCESSED  =  X0050D1 

OUTPUT  SAVED  IN  FILE  X0050D1.0UT 
DO  YOU  WANT  TO  MAKE  ANOTHER  RUN? 

ENTER  'YES'  OR  'NO', 
n 

*****  normal  termination  ***** 


Figure  15.  (Sheet  7  of  7) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.36.08 

**************** 

*  INPUT  DATA  * 

**************** 

1 . — HEADING 

'EXAMPLE  1  —  FIXED  END  BEAM  ANALYSIS 
'EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 


II. — BEAM  CROSS  SECTION  DATA 


X-COORDINATE 
START  STOP 

(IN)  (IN) 

.00  120.00 


MODULUS  OF 
ELASTICITY 
(PSI) 
3.00E+07 


<  - SECTION  PROPERTIES - > 

<  - START - >  < - STOP - > 

AREA  INERTIA  AREA  INERTIA 

(SQIN)  (IN**4)  (SQIN)  (IN**4) 

4.00  40.00  4.00  40.00 


. — NODE  SPACING  DATA 
X-COORDINATE 
START  STOP 

(FT)  (FT) 

-10.00  20.00 


MAXIMUM  NODE 
SPACING 
(FT) 

5.00 


IV. — LOAD  DATA 


IV. A. — CONCENTRATED  LOADS 


CONCENTRATED  LOADS 

X-COORD  X-LOAD  Y-LOAD  COUPLE 

(FT)  (K)  (K)  (K-FT) 

.00  40.00  .00  .00 


IV.  B.— DISTRIBUTED  LOADS 

LOAD  < - START - > 

DIRECT  X-COORD  LOAD 

(FT)  (K/FT) 

Y  .00  1.20 


X-COORD 

(FT) 

10.00 


LOAD 

(K/FT) 

1.20 


V. — FIXED  SUPPORT  DATA 
SUPPORT 
X-COORD 
(FT) 

.00 

10.00 


SPECIFIED  DISPLACEMENTS 
X-DISP.  Y-DISP.  ROTATION 

(FT)  (FT)  (RAD) 

FREE  .00  .00 

.00  .00  .00 


VI.  — LINEAR  SPRING  DATA 

NONE 

VII.  — NONLINEAR  SPRING  DATA 

NONE 


Chapter  5  Example  Solutions 
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■EXAMPLE  1  --  FIXED  END  BEAM  ANALYSIS 
'EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 
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Figure  17.  Axial  and  bending  stiffnesses  for  Example  1 


EXAMPLE  1  --  FIXED  END  BEAM  ANALYSIS 
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NONLINEAR  CONCENTRATED  SPRING 


■EXAMPLE  1  --  FIXED  END  BEAM  ANALYSIS 
•EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 
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Y-DISTRIBUTED  LOADS  (UNITS:  LB/IN.) 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.36.23 

************************ 

*  SUMMARY  OF  RESULTS  * 
************************ 

I. — HEADING 

'EXAMPLE  1  —  FIXED  END  BEAM  ANALYSIS 
'EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 


II.— MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM 

POSITIVE 

4.000E-02 

4.500E-02 

O.OOOE+00 

O.OOOE+00 

6.000E+03 

1.200E+05 


X-COORD 

(IN) 

.00 

60.00 

.00 

.00 

120.00 

.00 


MAXIMUM 
NEGATIVE 
O.OOOE+00 
O.OOOE+00 
O.OOOE+00 
-4.000E+04 
-6.000E+03 
-6. OOOE+04 


X-COORD 

(IN) 

.  00 
.00 
.00 
.00 
.00 
60.00 


III. — REACTIONS  AT  FIXED  SUPPORTS 
X-COORD 
(IN) 

.  00 

120.00 


X-REACTION 

(P) 

O.OOOE+00 

-4.000E+04 


Y-REACTION 

(P) 

-6.000E+03 

-6.000E+03 


MOM-REACTION 

(P-IN) 

-1.200E+05 

1.200E+05 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.37.34 

************************ 

*  SUMMARY  OF  RESULTS  * 
************************ 

I . — HEADING 

'EXAMPLE  lA  —  SAME  AS  EXAMPLE  1 

'EXCEPT  INCLUDE  EFFECTS  OF  AXIAL  LOAD  ON  BENDING 
'AXIAL  LOAD  =40  KIPS 


II. — MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM 
POSITIVE 
4 .OOOE-02 
4.555E-02 
O.OOOE+00 
O.OOOE+00 
6.000E+03 
1.210E+05 


X-COORD 

(IN) 

.00 

60.00 

.00 

.00 

120.00 

.00 


MAXIMUM 

X-COORD 

NEGATIVE 

(IN) 

O.OOOE+00 

.00 

O.OOOE+00 

.00 

O.OOOE+00 

.00 

-4.000E+04 

.00 

-6. OOOE+03 

.00 

-6.085E+04 

60.00 

III. — REACTIONS  AT  FIXED  SUPPORTS 

X-COORD  X-REACTION 

(IN)  (P) 

.00  O.OOOE+00 

120.00  -4.000E+04 


Y-REACTION 

(P) 

-6. OOOE+03 
-6. OOOE+03 


MOM-REACTION 
(P-IN) 
-1.210E+05 
1. 210E+05 


Figure  20.  Summary  of  results  for  initial  reruns  of  Example  1  (Sheet  1  of  3) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.38.08 

*★*★★★★★*★*★★★**★*★***★★ 

*  SUMMARY  OF  RESULTS  * 
★*★***★★★*★***★*****★**★ 

I.— HEADING 

'EXAMPLE  IB  —  INCREASE  AXIAL  LOAD  TO  2000  KIPS 


II . — MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM  X-COORD 


POSITIVE 
2 . OOOE+00 
1.138E-01 
O.OOOE+00 
0. OOOE+00 
6. OOOE+03 
2.373E+05 


(IN) 

.  00 
60.00 
.00 
.  00 

120.00 

.  00 


MAXIMUM 
NEGATIVE 
0 . OOOE+00 
0. OOOE+00 
0. OOOE+00 
-2 . OOOE+06 
-6. OOOE+03 
-1.703E+05 


X-COORD 
(IN) 
.00 
.00 
.00 
.00 
.  00 
60.00 


. — REACTIONS  AT  FIXED  SUPPORTS 

X-COORD  X-REACTION 

(IN)  (P) 

.00  O.OOOE+00 

120.00  -2.000E+06 


Y-REACTION  MOM-REACTION 

(P)  (P-IN) 

-6.000E+03  -2.373E+05 

-6.000E+03  2.373E+05 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.38.41 

*  SUMMARY  OF  RESULTS  * 

*★*★*★*★★*★★***★*★**★★★* 

I.— HEADING 

'EXAMPLE  1C  --  INCREASE  AXIAL  LOAD  TO  3000  KIPS 


II.— MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM  X-COORD 


POSITIVE 
3.000E+00 
5.041E-01 
0. OOOE+00 
O.OOOE+00 
6.000E+03 
8.818E+05 


(IN) 

.  00 
60.00 
.00 
.00 

120.00 

.00 


MAXIMUM 
NEGATIVE 
0. OOOE+00 
0. OOOE+00 
0. OOOE+00 
-3 . OOOE+06 
-6. OOOE+03 
-8. 103E+05 


X-COORD 

(IN) 

.00 

.00 

.00 

.00 

.00 

60.00 


. — REACTIONS  AT  FIXED  SUPPORTS 
X-COORD  > 

(IN) 

.  00 

120.00 


X-REACTION  Y-REACTION  MOM-REACTION 

(P)  (P)  (P-IN) 

O.OOOE+00  -6.000E+03  -8.818E+05 

-3.000E+06  -6.000E+03  8.818E+05 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-I994  TIME:  10.39.03 

*  SUMMARY  OF  RESULTS  * 
****★★****★**★********★* 

I.— HEADING 

'EXAMPLE  ID  —  INCREASE  AXIAL  LOAD  TO  3289  KIPS 
'THIS  AXIAL  LOAD  IS  SLIGHTLY  LESS  THAN  BUCKLING  LOAD 


II.  "MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

3 . 289E+00 

.  00 

O.OOOE+00 

.00 

1. 681E+02 

60.00 

O.OOOE+00 

.00 

0. OOOE+00 

.00 

O.OOOE+00 

.00 

O.OOOE+00 

.00 

-3.289E+06 

.00 

6.000E+03 

120.00 

-6. OOOE+03 

.00 

2.765E+08 

.00 

-2 .764E+08 

60.00 

III. — REACTIONS  AT  FIXED  SUPPORTS 
X-COORD 
(IN) 

.  00 

120.00 


X-REACTION 

(P) 

O.OOOE+00 

-3.289E+06 


Y-REACTION 

(P) 

-6. OOOE+03 
-6.000E+03 


MOM-REACTION 

(P-IN) 

-2.765E+08 

2.765E+08 


I 

I 


Figure  20.  (Sheet  3  of  3) 
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'EXAMPLE  1  —  FIXED  END  BEAM  ANALYSIS 
‘EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 
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Chapter  5  Example  Soultions 


LATERAL  DISPLACEMENTS  (UNITS:  IN) 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  10.39.12 

***************** 

*  INPUT  DATA  * 

★**★★★★*★***★★★★ 

I. — HEADING 

'EXAMPLE  IE  —  INCREASE  AXIAL  LOAD  TO  3290  KIPS 

'THIS  AXIAL  LOAD  IS  SLIGHTLY  GREATER  THAN  BUCKLING  LOAD 

II.  — BEAM  CROSS  SECTION  DATA 


< - SECTION  PROPERTIES- 


X-COORDINATE 
START  STOP 

(IN)  (IN) 

.00  120.00 


MODULUS  OF 
ELASTICITY 
(PSI) 
3.00E+07 


- START- 


AREA 

(SQIN) 

4.00 


INERTIA 

(IN**4) 

40.00 


< - STOP - > 

AREA  INERTIA 

(SQIN)  (IN**4) 

4.00  40.00 


III.  — NODE  SPACING  DATA 

X-COORDINATE  MA 

START  STOP 

(FT)  (FT) 

-10.00  20.00 

IV. — LOAD  DATA 

IV . A . — CONCENTRATED  LOADS 


MAXIMUM  NODE 
SPACING 
(FT) 

5.00 


CONCENTRATED  LOADS 


COORD 

X-LOAD 

Y-LOAD 

(FT) 

(K) 

(K) 

.00 

40.00 

.00 

.00 

1960.00 

.00 

.00 

1000.00 

.00 

.00 

289.00 

.00 

.00 

1.00 

.00 

IV. B. — DISTRIBUTED  LOADS 


LOAD 

DIRECT 


- START - > 

-COORD  LOAD 

(FT)  (K/FT) 

.00  1.20 


COUPLE 

(K-FT) 

.00 
.  00 
.00 
.00 
.  00 

< - STOP- 

X-COORD 

(FT) 

10.00 


- > 

LOAD 

(K/FT) 

1.20 


V.  — FIXED  SUPPORT  DATA 

SUPPORT 

X-COORD  X-E 

(FT)  (F 

.  00  FF 

10.00 

VI.  — LINEAR  SPRING  DATA 

NONE 

VII.  — NONLINEAR  SPRING  DATA 

NONE 


SPECIFIED  DISPLACEMENTS 
-DISP.  y-DISP.  ROTATION 

[FT)  (FT)  (RAD) 

■■REE  .00  .00 

.00  .00  .00 
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PROGRAM  CBEAMC 

-  ANALYSIS 

OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 

DATE:  21-MAR-1994 

TIME:  10.39.30 

************************ 

*  SUMMARY  OF  RESULTS  * 

************************ 

I .  — HEADING 

'EXAMPLE  IE  - 

-  INCREASE 

AXIAL  LOAD  TO  3290 

KIPS 

'THIS  AXIAL  LOAD  IS  SLIGHTLY  GREATER  THAN 

BUCKLING  LOAD 

II. — MAXIMA 

MAXIMUM  X 

-COORD  MAXIMUM 

X-COORD 

POSITIVE 

(IN)  NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN) 

3-290E+00 

.00  O.OOOE+00 

.00 

LATERAL  DISPLACEMENT  (IN) 

O.OOOE+00 

.00  -1.106E+03 

60.00 

ROTATION  (RAD) 

O.OOOE+00 

.00  O.OOOE+00 

.00 

AXIAL  FORCE  (P) 

O.OOOE+00 

.00  -3.290E+06 

.00 

SHEAR  (P) 

6.000E+03 

120.00  -6.000E+03 

.00 

BENDING  MOMENT 

(P-IN) 

1.820E+09 

60.00  -1.820E+09 

.00 

III. — REACTIONS 

AT  FIXED  SUPPORTS 

X-COORD 

X-REACTION 

Y-REACTION  MOM-REACTION 

(IN) 

(P) 

(P)  (P- 

IN) 

.00 

O.OOOE+00 

-6.000E+03  1.820E+09 

120.00 

-3.290E+06 

-6.000E+03  -1.820E+09 

Figure  23.  (Concluded) 


Example  2:  Beam  on  Uniform  Elastic  Foundation 


Accuracy  of  results  produced  by  the  finite  element  model  used  to  represent 
the  system  shown  in  Figure  24  depends  on  node  spacing.  Solutions  were 
obtained  for  two  node  spacings  as  indicated  in  the  data  file  shown  in 
Figure  25. 


The  echoprint  of  input  data,  summary  of  results,  and  the  (trivial)  complete 
results  table  for  the  first  case  Oines  1000  through  1140  of  the  input  file)  are 
shown  in  Figure  26.  Maximum  lateral  deflection  and  bending  moment  for  the 
model  with  three  nodes  and  two  elements  differ  only  slightly  from  the  “exact” 
solution  (0.19097  in.  and  -1.52025E+5  lb-in.,  respectively). 

Increasing  the  model  to  11  nodes  and  10  elements  produces  results.  Fig¬ 
ure  27,  which  are  exact  to  four  significant  figures. 

It  should  be  emphasized  that  node  spacing  may  have  a  more  significant 
effect  on  the  accuracy  of  the  solution  than  is  indicated  by  this  problem.  This 
is  demonstrated  in  Example  3. 


45 


1.2  k/ft 


?  M  n  Ml' 

k  =100  psi 

y 


Cross  Section:  E  =  3010^ 
A  =  4  sq  in. 
1  =  40  in.^ 


Figure  24.  Beam  on  elastic  support  for  Example  2 
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1000 

'EXAMPLE  2  ~  SIMPLE 

BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 

1010 

'ILLUSTRATE  EFFECT 

OF  NODE  SPACING  ON  SOLUTION 

1020 

'NODES  AT  5  FT. 

1030 

BEAM 

1040 

0  120  30. E6  4  40 

1050 

NODES  F 

1060 

0  10  5 

1070 

LOADS  F  K 

1080 

D  Y  0  1.2  10 

1090 

FIXED  F 

1100 

0  FREE  0  FREE 

1110 

10  0  0  FREE 

1120 

LINEAR 

1130 

D  Y  0  120  100  0  0 

1140 

FINISH  RERUN 

1150 

'EXAMPLE  2 A  —  SAME 

AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 

1160 

NODES  N  F 

1170 

0  10  1.25 

1180 

FINISH 

Figure  25.  Input  file  for  Example  2 


Chapter  5  Example  Solutions 


47 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.04.36 

iric‘kicit1tie'k1e'kifk*-k-kic 

*  INPUT  DATA  * 

•kitifkifkifkiticicitic-kic* 

I .  — HEADING 

'EXAMPLE  2  -  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
'ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
'NODES  AT  5  FT. 

II.  — BEAM  CROSS  SECTION  DATA 

< - SECTION  PROPERTIES - > 

X-COORDINATE  MODULUS  OF  < - START - >  < - STOP - > 

START  STOP  ELASTICITY  AREA  INERTIA  AREA  INERTIA 

(IN)  (IN)  (PSI)  (SQIN)  (IN**4)  (SQIN)  (IN**4) 

.00  120.00  3.00E+07  4.00  40.00  4.00  40.00 

III.  — NODE  SPACING  DATA 

X-COORDINATE  MAXIMUM  NODE 

START  STOP  SPACING 

(FT)  (FT)  (FT) 

.00  10.00  5.00 

IV. — LOAD  DATA 

4. A. — CONCENTRATED  LOADS 
NONE 

IV. B. — DISTRIBUTED  LOADS 

LOAD  < - START - >  < - STOP - > 

DIRECT  X-COORD  LOAD  X-COORD  LOAD 

(FT)  (K/FT)  (FT)  (K/FT) 

Y  .00  1.20  10.00  1.20 

V. — FIXED  SUPPORT  DATA 

SUPPORT  SPECIFIED  DISPLACEMENTS 

X-COORD  X-DISP.  Y-DISP.  ROTATION 

(FT)  (FT)  (FT)  (RAD) 

•00  FREE  .00  FREE 

10.00  .00  .00  FREE 

VI. — LINEAR  SPRING  DATA 

VI. A. — CONCENTRATED  LINEAR  SPRINGS 
NONE 

VI.  B.— DISTRIBUTED  LINEAR  SPRINGS 

SPRING  X-COORD  SPRING  STIFFNESS  COEFFICIENTS 

DIRECTION  START  STOP  ABC 

(IN)  (IN)  (P/I/I)  (P/I/I) 

Y  .00  120,00  100.00  .00  .00 

VII.  — NONLINEAR  SPRING  DATA 

NONE 


Figure  26.  Echoprint  of  input  and  results  for  Example  2A  (Continued) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.04.53 

it****lt1fk*'kit*-kifitic*it*'k**** 

*  SUMMARY  OF  RESULTS  * 

it'kit'kic'kitit'kifkiticitlticitit'klticit'kls 

I ,  —HEADING 

'EXAMPLE  2  -  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
'ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
'NODES  AT  5  FT. 

II .  —MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  fP) 

BENDING  MOMENT  (P-IN) 


MAXIMUM 

POSITIVE 

O.OOOE+00 

1.912E-01 

5.116E-03 

O.OOOE+00 

5.273E+03 

O.OOOE+00 


X-COORD 

(IN) 

.00 

60.00 

.00 

.00 

120.00 


MAXIMUM 

NEGATIVE 

O.OOOE+00 

O.OOOE+00 

-5.116E-03 

O.OOOE+OO 

-5.273E+03 

-1.522E+05 


III . —REACTIONS  AT  FIXED  SUPPORTS 

X-COORD  X-REACTION  Y-REACTION 

(IN)  (P)  (P) 

.00  O.OOOE+OO  -5.273E+03 

120.00  O.OOOE+OO  -5.273E+03 


X-COORD 

(IN) 

.00 

.00 

120.00 

.00 

.00 

60.00 


MOM-REACTION 

(P-IN) 

O.OOOE+OO 

O.OOOE+OO 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.04.53 

*  COMPLETE  RESULTS  * 

iiltitititicit'itititifk'k'kifkifk'kititit 

I .  —HEADING 

'EXAMPLE  2  -  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
'ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
'NODES  AT  5  FT. 

II. — DISPLACEMENTS  AND  INTERNAL  FORCES 

< - DISPLACEMENTS - >  < - INTERNAL  FORCES- 


X-COORD  AXIAL  LATERAL  ROTATION 

(IN)  (IN)  (IN)  (RAD) 

.00  O.OOOE+OO  O.OOOE+OO  5.116E-03 

60.00  O.OOOE+OO  1.912E-01  1.542E-19 

120.00  O.OOOE+OO  O.OOOE+OO  -5.116E-03 

III. — FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 


AXIAL  SHEAR  MOMENT 

(P)  (P)  (P-IN) 

O.OOOE+OO  -5.273E+03  O.OOOE+OO 
O.OOOE+OO  O.OOOE+OO  -1.522E+05 
O.OOOE+OO  5.273E+03  O.OOOE+OO 


X-COORD 

(IN) 

.00 

60.00 

120.00 


AXIAL 

(P/IN) 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 


LATERAL 

(P/IN) 

O.OOOE+OO 

-1.912E+01 

O.OOOE+OO 


Figure  26.  (Concluded) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.05.35 

**************** 

*  INPUT  DATA  * 

**************** 

I .  — HEADING 

'EXAMPLE  2A  —  SAME  AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 

II.  — BEAM  CROSS  SECTION  DATA 


X-COORDINATE  MODULUS  OF 

START  STOP  ELASTICITY 

(IN)  (IN)  (PSI) 

.00  120.00  3.00E+07 

III. — NODE  SPACING  DATA 

X-COORDINATE 
START  STOP 

(FT)  (FT) 

.00  10.00 

IV.  — LOAD  DATA 
4. A. — CONCENTRATED  LOADS 

NONE 

IV. B. — DISTRIBUTED  LOADS 

LOAD  < - START - > 

DIRECT  X-COORD  LOAD 

(FT)  (K/FT) 

Y  .00  1.20 

V.  — FIXED  SUPPORT  DATA 

SUPPORT 
X-COORD 
(FT) 

.00 

10.00 

VI. — LINEAR  SPRING  DATA 
VI. A.— CONCENTRATED  LINEAR  SPRINGS 

NONE 

VI. B. — DISTRIBUTED  LINEAR  SPRINGS 

SPRING  X-COORD  SPRING  STIFFNESS  COEFFICIENTS 

DIRECTION  START  STOP  ABC 

(IN)  (IN)  (P/I/I)  (P/I/I) 

Y  .00  120.00  100.00  .00  .00 

VII.  — NONLINEAR  SPRING  DATA 

NONE 


Figure  27.  Echoprint  of  input  and  results  for  Example  2B  (Sheet  1  of  3) 


SPECIFIED  DISPLACEMENTS 
X-DISP.  Y-DISP.  ROTATION 

(FT)  (FT)  (RAD) 

FREE  .00  FREE 

.00  .00  FREE 


< - STOP - > 

X-COORD  LOAD 

(FT)  (K/FT) 

10.00  1.20 


MAXIMUM  NODE 
SPACING 
(FT) 

1.25 


<  - SECTION  PROPERTIES - > 

<  - START - >  < - STOP - > 

AREA  INERTIA  AREA  INERTIA 

(SQIN)  (IN**4)  (SQIN)  (IN**4) 

4.00  40.00  4.00  40.00 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLJNE^ 
DATE:  21-MAR-1994 

************************ 

*  SUMMARY  OF  RESULTS  * 
************************ 


SUPPORTS 

16.05.43 


-  SAME  AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 


II .  —MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

O.OOOE+OO 

.00 

O.OOOE+OO 

.00 

1.910E-01 

60.00 

O.OOOE+OO 

,  00 

5.109E-03 

.00 

-5.109E-03 

120.00 

O.OOOE+OO 

.00 

O.OOOE+OO 

.00 

5.266E+03 

120.00 

-5,266E+03 

,  00 

O.OOOE+OO 

,00 

-1.520E+05 

60.00 

III.— reactions  AT  FIXED  SUPPORTS 


X-COQRD 

(IN) 

.00 

120.00 


X-REACTION 

(P) 

O.OOOE+00 

O.OOOE+00 


y-REACTION  MOM-REACTION 
(P)  (P-IN) 

-5.266E+03  O.OOOE+OO 

-5.266E+03  O.OOOE+OO 


Figure  27.  (Sheet  2  of  3) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.05.43 


********************** 


I.  "HEADING 
'EXAMPLE  2A 


*  COMPLETE  RESULTS  * 

****************«*«*^* 


SAME  AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 


II. — DISPLACEMENTS  AND  INTERNAL  FORCES 
< - DISPLACEMENTS - > 


X-COORD 

AXIAL 

(IN) 

(IN) 

.00 

O.OOOE+00 

15.00 

0. OOOE+00 

30.00 

O.OOOE+00 

45.00 

O.OOOE+00 

60.00 

O.OOOE+00 

75.00 

O.OOOE+00 

90.00 

O.OOOE+00 

105.00 

O.OOOE+00 

120.00 

O.OOOE+00 

LATERAL 

ROTATION 

(IN) 

(RAD) 

0. OOOE+00 

5.109E-03 

7.434E-02 

4.661E-03 

1.362E-01 

3.495E-03 

1.769E-01 

1.862E-03 

1.910E-01 

-4.597E-18 

1.769E-01 

-1.862E-03 

1.362E-01 

-3.495E-03 

7.434E-02 

-4.661E-03 

0. OOOE+00 

-5.109E-03 

< - INTERNAL  FORCES 


AXIAL 

(P) 

0. OOOE+00 
0. OOOE+00 
O.OOOE+00 
0. OOOE+00 
O.OOOE+00 
0. OOOE+00 
O.OOOE+00 
O.OOOE+00 
O.OOOE+00 


SHEAR 

(P) 

-5.266E+03 

-3.823E+03 

-2.483E+03 

-1.221E+03 

O.OOOE+00 

1.221E+03 

2.483E+03 

3.823E+03 

5.266E+03 


MOMENT 

(P-IN) 

0. OOOE+00 
-6.802E+04 
-1. 152E+05 
-1.429E+05 
-1.520E+05 
-1.429E+05 
-1.152E+05 
-6.802E+04 
O.OOOE+00 


III." FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 


DISTRIBUTED  SPRING  FORCES 


X-COORD 

AXIAL 

LATERAL 

(IN) 

(P/IN) 

(P/IN) 

.  00 

0. OOOE+00 

O.OOOE+OO 

15.00 

O.OOOE+00 

-7.434E+00 

30.00 

0. OOOE+00 

-1.362E+01 

45.00 

O.OOOE+00 

-1.769E+01 

60.00 

0. OOOE+00 

-1.910E+01 

75.00 

O.OOOE+00 

-1.769E+01 

90.00 

O.OOOE+OO 

-1.362E+01 

105.00 

0. OOOE+OO 

-7.434E+00 

120.00 

O.OOOE+OO 

O.OOOE+OO 

Figure  27.  (Sheet  3  of  3) 
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Example  3:  Pile  Head  Stiffness  Matrix 


The  program  was  used  to 
generate  the  pile  head  stiffness 
matrix  for  the  system  shown  in 
Figure  28.  By  definition,  coeffi¬ 
cients  of  the  pile  head  stiffness 
matrix  are  the  reactions  gener¬ 
ated  by  imposing  unit  values  of 
deflection  components  at  the  pile 
head. 

The  problem  was  run  using 
the  input  data  file  shown  in  Fig¬ 
ure  29.  Because  axial  load 
effects  on  bending  were 
excluded,  axial  and  lateral  coeffi¬ 
cients  for  the  stiffness  matrix  are 
obtained  simultaneously  by 
imposing  unit  axial  (x)  and 
lateral  (y)  deflections  at  the  pile 
head  with  rotation  at  that  point 
equal  to  zero.  Two  solutions  are 
presented  for  this  case;  for 
nodes  spaced  at  5  ft  and  nodes 
spaced  at  1  ft.  Moment  coeffi¬ 
cients  for  the  pile  head  stiffness 
matrix  were  obtained  for  a  1-ft 
node  spacing  by  specifying  zero 
axial  and  lateral  displacements 
and  a  unit  rotation  at  the  pile 
head. 


Results  for  each  case  are 
given  in  Figures  30  and  31. 

Values  of  stiffness  coefficients  obtained  for  the  5-ft  node  spacing  are  signifi¬ 
cantly  in  error.  Decreasing  the  node  spacing  to  1  ft  produced  the  pile  head 
stiffness  matrix  shown  as  follows.  Values  shown  in  parentheses  were  obtained 
by  another  approximate  method  reported  by  Reese  and  Matlock  (1956). 


Figure  28.  Pile/soil  system  for  Example  3 


1 

- 

“ 

1 

3.602E6 

0 

0 

2.Q11E5 

5A2AE6 

u 

'Py 

* 

(2.014E5) 

{5.356E6) 

4 

V  • 

M, 

SYM 

2.148E8 

0 

V.  J 

(2.137ES). 

Reducing  the  node  spacing  to  0.5  ft  had  little  effect  on  the  results. 
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1000  'EXAMPLE  3  —  PILE  STIFFNESS  MATRIX 
1010  'AXIAL  FORCE  EFFECTS  OF  BENDING  EXCLUDED 
1020  'NODES  AT  5  FT. 

1030  'AXIAL  AND  LATERAL  TERMS 
1040  BEAM 

1050  0  1200  4.3E6  100  833.33 

1060  NODES  F 

1070  0  100  5 

1080  FIXED  I 

1090  0110 

1100  LINEAR 

1110  D  Y  0  1200  0  10  2 

1120  D  X  0  1200  0  10  2 

1130  FINISH  RERUN 

1140  'EXAMPLE  3  —  CONTINUED 

1150  'AXIAL  AND  LATERAL  TERMS 

1160  'NODES  AT  1  FT. 

1170  NODES  NEW  F 
1180  0  100  1 
1190  FINISH  RERUN 
1200  'EXAMPLE  3  —  CONTINUED 
1210  'MOMENT  TERMS 
1220  FIXED  NEW 
1230  0001 
1240  FINISH 


Figure  29.  Input  file  for  Example  3 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.22.10 

ic'kltit’kititltititit'k'kic'k'k 

*  INPUT  DATA  * 

ifkieit-kifk'kltifkit'kic'kic 

I .  —HEADING 

'EXAMPLE  3  —  PILE  STIFFNESS  MATRIX 
'AXIAL  FORCE  EFFECTS  OF  BENDING  EXCLUDED 
'NODES  AT  5  FT. 

'XIAL  AND  LATERAL  TERMS 

II.  — BEAM  CROSS  SECTION  DATA 

< - SECTION  PROPERTIES - > 

X-COORDINATE  MODULUS  OF  < - START - >  < - STOP - > 

START  STOP  ELASTICITY  AREA  INERTIA  AREA  INERTIA 

(IN)  (IN)  (PSI)  (SQIN)  (IN**4)  (SQIN)  (IN**4) 

.00  1200.00  4.30E+06  100.00  833.33  100.00  833.33 

III.  — NODE  SPACING  DATA 

X-COORDINATE  MAXIMUM  NODE 

START  STOP  SPACING 

(FT)  (FT)  (FT) 

.00  100.00  5.00 

IV.  — LOAD  DATA 

NONE 

V.  — FIXED  SUPPORT  DATA 

SUPPORT 
X-COORD 
(IN) 

.00 

VI. — LINEAR  SPRING  DATA 
VI. A. — CONCENTRATED  LINEAR  SPRINGS 

NONE 

VI. B. — DISTRIBUTED  LINEAR  SPRINGS 
SPRING  X-COORD 

DIRECTION  START  STOP 

(IN)  (IN) 

Y  .00  1200.00 

X  .00  1200.00 

VII.  — NONLINEAR  SPRING  DATA 

NONE 


Figure  30.  Output  for  Example  3,  model  with  5-ft  node  spacing  (Continued) 


SPRING  STIFFNESS  COEFFICIENTS 
ABC 
(P/I/I)  (P/I/I) 

.00  10.00  2.00 

.00  10.00  2.00 


SPECIFIED  DISPLACEMENTS 
X-DISP.  Y-DISP.  ROTATION 

(IN)  (IN)  (RAD) 

1.00  1.00  .00 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.22.23 

*  SUMMARY  OF  RESULTS  * 

'kiticicititit'k'kitififkititititititifkititit 

I. — HEADING 

'EXAMPLE  3  —  PILE  STIFFNESS  MATRIX 
'AXIAL  FORCE  EFFECTS  OF  BENDING  EXCLUDED 
'NODES  AT  5  FT. 

'XIAL  AND  LATERAL  TERMS 

II .  —MAXIMA 

MAXIMUM  X-COORD  MAXIMUM  X-COORD 

POSITIVE  (IN)  NEGATIVE  (IN) 

AXIAL  DISPLACEMENT  (IN)  :  l.OOOE+00  .00  -8.297E-06  300.00 

LATERAL  DISPLACEMENT  (IN):  l.OOOE+00  .00  -l,547E-04  180.00 

ROTATION  (RAD)  :  3.749E-05  120.00  -6.680E-03  60.00 

AXIAL  FORCE  (P)  :  2.234E+02  300.00  -3.892E+06  .00 

SHEAR  (P)  :  3.537E+05  .00  -5.108E+04  60.00 

BENDING  MOMENT  (P-IN)  :  1.690E+06  60.00  -7,331E+06  .00 

III.  — REACTIONS  AT  FIXED  SUPPORTS 

X-COORD  X-REACTION  Y-REACTION  MOM-REACTION 

(IN)  (P)  (P)  (P-IN) 

.00  3.892E+06  3.537E+05  7.331E+06 

PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.22.23 

********************** 

*  COMPLETE  RESULTS  * 

********************** 

I.  — HEADING 

'EXAMPLE  3  —  PILE  STIFFNESS  MATRIX 
'AXIAL  FORCE  EFFECTS  OF  BENDING  EXCLUDED 
'NODES  AT  5  FT. 

'AXIAL  AND  LATERAL  TERMS 

II. — DISPLACEMENTS  AND  INTERNAL  FORCES 


X-COORD 

AXIAL 

DISPLACEMENT 

LATERAL 

ROTATION 

< - 

AXIAL 

JTERNAL  FORC 
SHEAR 

:es - > 

MOMENT 

(IN) 

(IN) 

(IN) 

(RAD) 

(P) 

(P) 

(P-IN) 

.00 

l.OOOE+00 

l.OOOE+00 

0. OOOE+00 

-3.892E+06 

3.537E+05 

-7.331E+06 

60.00 

4.944E-01 

4.959E-02  -6.680E-03 

-3.176E+06 

-5.108E+04 

1.690E+06 

120.00 

1.580E-01 

7.693E-04 

3.749E-05 

-1.596E+06 

4.434E+03 

-5.210E+04 

III.— FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 
X-COORD  AXIAL  LATERAL 

(IN)  (P/IN)  (P/IN) 

.00  O.OOOE+00  O.OOOE+00 

60.00  -1.780E+04  -1.785E+03 

120.00  -2.276E+04  -1.108E+02 


Figure  30.  (Concluded) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.23.32 

************************ 

*  SUMMARY  OF  RESULTS  * 

************************ 

I .  —HEADING 

'EXAMPLE  3  —  CONTINUED 
'AXIAL  AND  LATERAL  TERMS 
'NODES  AT  1  FT. 

II .  —MAXIMA 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN) 

l.OOOE+00 

.00 

0. OOOE+00 

.00 

LATERAL  DISPLACEMENT  (IN) 

l.OOOE+00 

.00 

-1.716E-02 

96.00 

ROTATION  (RAD) 

6.306E-04 

108.00 

-2.012E-02 

24.00 

AXIAL  FORCE  (P) 

O.OOOE+00 

,00 

-3.602E+06 

.00 

SHEAR  (P) 

2.077E+05 

.00 

-4.836E+04 

72.00 

BENDING  MOMENT  (P-IN) 
III.— REACTIONS  AT  FIXED  £ 

1.709E+06 

SUPPORTS 

60,00 

-5.424E+06 

.00 

X-COORD 

X-REACTION 

Y-REACTION  MOM-REACTION 

(IN) 

(P) 

(P)  (P- 

IN) 

.00 

3.602E+06 

2.077E+05  5.424E+06 

PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME:  16.23.32 

Ititit’kitititieilititititititit'kititit'kit 

*  COMPLETE  RESULTS  * 

********************** 

I . — HEADING 

'EXAMPLE  3  —  CONTINUED 
'AXIAL  AND  LATERAL  TERMS 
'NODES  AT  1  FT. 


II . —DISPLACEMENTS  AND  INTERNAL  FORCES 


< - DIS  PLACEMENTS - > 

< - INTERNAL  FOR( 

:es - > 

:-COORD 

AXIAL 

LATERAL 

ROTATION 

AXIAL 

SHEAR 

MOMENT 

(IN) 

(IN) 

(IN) 

(RAD) 

(P) 

(P) 

(P-IN) 

.00 

l.OOOE+00 

l.OOOE+00 

O.OOOE+00 

-3.602E+06 

2.077E+05 

-5.424E+06 

12.00 

8.996E-01 

9.076E-01 

-1.402E-02 

-3.594E+06 

1.994E+05 

-2.966E+06 

24.00 

7,997E-01 

6.953E-01 

-2.012E-02 

-3.557E+06 

1.655E+05 

-7.442E+05 

36.00 

7.015E-01 

4.511E-01 

-1.967E-02 

-3.474E+06 

1.029E+05 

8.879E+05 

48.00 

6. 063E-01 

2.398E-01 

-1.511E-02 

-3.334E+06 

3.139E+04 

1.689E+06 

60.00 

5.159E-01 

9.371E-02 

-9.233E-03 

-3.136E+06 

-2.371E+04 

1.709E+06 

72.00 

4.318E-01 

1.465E-02 

-4.205E-03 

-2.888E+06 

-4.836E+04 

1.245E+06 

84.00 

3 .551E-01 

-1.476E-02 

-1.026E-03 

-2.600E+06 

-4.545E+04 

6.602E+05 

96.00 

2.869E-01 

-1.716E-02 

3.740E-04 

-2.288E+06 

-2.845E+04 

2.104E+05 

108.00 

2.275E-01 

-1.036E-02 

6.306E-04 

-1.967E+06 

-1. 106E+04 

-2.207E+04 

120.00 

1.770E-01 

-3.864E-03 

4.187E-04 

-1.652E+06' 

-4 . 502E+02 

-8.321E+04 

Chapter  5  Example  Solutions 


III . — FORCES 

IN  DISTRIBUTED  LINEAR  SPRINGS 

DISTRIBUTED  SPRING  FORCES 

X-COORD 

AXIAL 

LATERAL 

(IN) 

(P/IN) 

(P/IN) 

.  00 

0. OOOE+00 

O.OOOE+00 

12 . 00 

-l,295E+03 

-1. 307E+03 

24 . 00 

-4.606E+03 

-4 . 005E+03 

36.00 

-9.091E+03 

-5.847E+03 

48.00 

-1.397E+04 

-5.525E+03 

60.00 

-1.857E+04 

-3 . 373E+03 

72 . 00 

-2.238E+04 

-7.596E+02 

84 . 00 

-2.506E+04 

1, 042E+03 

96.00 

-2.644E+04 

1.582E+03 

108.00 

-2.654E+04 

1.208E+03 

120.00 

-2.549E+04 

5.564E+02 

Figure  31 .  (Sheet  2  of  4) 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE;  21-MAR-1994  TIME:  16.23.49 

it*********************** 

*  SUMMARY  OF  RESULTS  * 

★  ★************<t********* 

I .  —HEADING 

'EXAMPLE  3  —  CONTINUED 
'MOMENT  TERMS 


II. — MAXIMA 


MAXIMUM  X- 

-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN) 

O.OOOE+00 

.00 

O.OOOE+00 

.00 

LATERAL  DISPLACEMENT  (IN) 

1.021E+01 

24.00 

-2.083E-01 

108.00 

ROTATION  (RAD) 

l.OOOE+00 

.00 

-2.501E-01 

48.00 

AXIAL  FORCE  (P) 

O.OOOE+00 

.00 

O.OOOE+00 

.00 

SHEAR  (P) 

5.424E+06 

.00 

-7.242E+05 

84.00 

BENDING  MOMENT  (P-IN) 

2.380E+07 

72.00 

-2.148E+08 

.00 

III. — REACTIONS  AT  FIXED  SUPPORTS 

X-COORD 

X-REACTION 

Y-REACTION  MOM-REACTION 

(IN) 

(P) 

(P) 

(P- 

IN) 

.00 

O.OOOE+00 

5.424E+06  2.148E+08 

PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  21-MAR-1994  TIME;  16.23.49 

1t‘k-k‘k‘kir1t**it‘k‘krk-k1tit‘k‘k'kit‘k-k 

*  COMPLETE  RESULTS  * 

Itifk’kit'kifkititleitifk-kifkit'kiticit 

I . — HEADING 

'EXAMPLE  3  —  CONTINUED 
'MOMENT  TERMS 


II.— 

•DISPLACEMENTS 

AND  INTERNAL  FORCES 

X-COORD 

AXIAL 

LATERAL 

ROTATION 

AXIAL 

SHEAR 

MOMENT 

(IN) 

(IN) 

(IN) 

(RAD) 

(P) 

(P) 

(P-IN) 

.00 

0. OOOE+00 

O.OOOE+OO 

l.OOOE+00 

O.OOOE+OO 

5.424E+06 

-2.148E+08 

12.00 

0. OOOE+00 

8.121E+00 

3.897E-01 

O.OOOE+OO 

5.372E+06 

-1.498E+08 

24.00 

0. OOOE+00 

1.021E+01 

-5.985E-03 

O.OOOE+OO 

4.950E+06 

-8.732E+07 

36.00 

0. OOOE+00 

8.767E+00 

-2.050E-01 

O.OOOE+OO 

3.873E+06 

-3.372E+07 

48.00 

0. OOOE+00 

5.911E+00 

-2.501E-01 

O.OOOE+OO 

2. 307E+06 

3.630E+06 

60.00 

O.OOOE+00 

3.136E+00 

-2.024E-01 

O.OOOE+OO 

7.573E+05 

2. 173E+07 

72.00 

O.OOOE+00 

1.179E+00 

-1.226E-01 

O.OOOE+OO 

-3.067E+05 

2.380E+07 

84 . 00 

0 . COOE+00 

1.494E-01 

-5.290E-02 

O.OOOE+OO 

-7.242E+05 

1.700E+07 

96.00 

O.OOOE+00  - 

■2.027E-01 

-1.062E-02 

O.OOOE+OO 

-6.495E+05 

8.407E+06 

108.00 

O.OOOE+OO  - 

■2.083E-01 

6.225E-03 

O.OOOE+OO 

-3.709E+05 

2.220E+06 

120.00 

O.OOOE+00  - 

■1.132E-01 

8.053E-03 

O.OOOE+OO 

-1.192E+05 

-6.218E+05 

Figure  31.  (Sheet  3  of  4) 
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III. — FORCES 

IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 

X-COORD 

AXIAL 

LATERAL 

(IN) 

(P/IN) 

(P/IN) 

.00 

0. OOOE+00 

O.OOOE+00 

12 . 00 

O.OOOE+00 

-1. 169E+04 

24 . 00 

0. OOOE+00 

-5.883E+04 

36.00 

O.OOOE+00 

-1.136E+05 

48 . 00 

O.OOOE+00 

-1. 362E+05 

60.00 

0. OOOE+00 

-1.129E+05 

72.00 

O.OOOE+00 

-6.113E+04 

84.00 

O.OOOE+00 

-1. 054E+04 

96.00 

O.OOOE+00 

1.868E+04 

108.00 

O.OOOE+00 

2.430E+04 

120.00 

O.OOOE+00 

1. 630E+04 

Figure  31 .  (Sheet  4  of  4) 


Example  4:  Anchored  Retaining  Wall 

A  1-ft  strip  of  the  multiple  anchored  retaining  wall  described  in  Figure  32 
was  analyzed  as  a  beam-column  in  which  the  soil  and  anchors  were  repre¬ 
sented  by  nonlinear  springs. 

Procedures  described  by  Haliburton  (1971)  were  used  to  obtain  the  force- 
deformation  curves  shown  in  Figure  33  for  the  soil-resisting  lateral  (y)  dis¬ 
placements  of  the  system.  To  the  writer’s  knowledge,  no  procedure  has  been 
established  for  representing  wall-soil  friction  interaction.  The  force- 
displacement  curves  shown  in  Figure  34  were  arbitrarily  selected  for  this 
effect. 

The  anchors  were  represented  by  concentrated  nonlinear  springs  with  char¬ 
acteristics  shown  in  Figure  35. 

The  input  data  file  used  for  solution  with  the  program  is  shown  in  Fig¬ 
ure  36.  Output  from  the  program  is  given  in  Figure  37.  Anchor  forces  are 
given  by  Section  III  of  the  Summary  of  Results.  Soil  pressures  are  given  by 
the  tabulation  of  forces  in  distributed  nonlinear  springs  in  Section  III  of  the 
Complete  Results.  Plots  of  the  input  data  are  shown  in  Figures  38  through  40. 
Plots  of  the  results  are  shown  in  Figures  41  through  43. 
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Beam:  PZ38  Sheet  Pile 
E  =  29-10®  psi 
A  =  11.2  sq  ia/ft  of  wall 
1  =  280.0  ia^  /ft  of  wall 


Soil  1 

Soil  2 

Unit  Weight 

100  pcf 

62.5  pcf  (sub; 

Internal  Friction 

o 

O 

CO 

0 

Cohesion 

0 

1000  psf 

At  Rest,  ko 

0.5 

1 

Soil  Modulus 

8  pci 

50  pci 

Interact.  Dist. 

6ft 

10ft 

Anchor  1 

Anchor  2 

E 

29  -10®  psi 

29  -10®  psi 

A  (sq  in. /ft) 

0.16 

0.2 

Eff.  Length  (ft) 

25 

30 

1 

Figure  32.  Multiple  anchored  retaining  wall  for  Example  4 
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Figure  33.  Force-deformation  curves  for  lateral  (y)  soil  resistance 
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Leftside  at  x  =  26' 


Leftside  at  x  =  46 


Figure  34.  Force-deformation  curves  for  axial  (x)  soil  resistance 
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Force  Parallel 
to  y-Axis  (kips) 


Anchor  at  x  =  6' 


Force  at  70  ° 
with  x-Axis  (kips) 


Anchor  at  x  =  17' 


Figure  35.  Force-deformation  curves  for  anchors 
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1000  'EXAMPLE  4  —  MULTIPLE  ANCHORED  RETAINING  WALL 
1010  'EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
1020  'TENSION-ONLY  ANCHORS 
1030  BEAM 

1040  (PZ  38  SHEET  PILE) 

1050  0  552  29. E6  11.2  280.8 
1060  NODES  F 
1070  0  50  1 
1080  NONLINEAR  F  P 

1090  (LATERAL  SPRINGS  FOR  RIGHTSIDE  SOIL) 

1100  D  Y  0  2  1  1 

1110  -0.008  -133  0.125  -1200 

1120  C  16  2  1  1 

1130  -0.008  -747  0.125  -6720 

1140  C  17  2  1  1 

1150  -0.231  -303  0.231  -4303 

1160  E  46  2  1  1 

1170  -0.231  -2115  0.231  -6115 

1180  (LATERAL  SPRINGS  FOR  LEFTSIDE  SOIL) 

1190  D  Y  26  2  1  1 

1200  -0.231  2000  0.231  0 

1210  E  46  2  1  1 

1220  -0.231  3250  0.145  0 

1230  (AXIAL  SPRINGS  FOR  RIGHTSIDE  SOIL) 

1240  D  X  0  2  1  1 

1250  0  0  0.1  -61 

1260  C  16  2  1  1 

1270  0  0  0.1  -342 

1280  C  17  2  1  1 

1290  -0.1  500  0.1  -500 

1300  E  46  2  1  1 

1310  -0.1  1000  0.1  -1000 

1320  (AXIAL  SPRINGS  FOR  LEFTSIDE  SOIL) 

1330  D  X  26  2  1  1 
1340  -0.1  500  0.1  -500 
1350  E  46  2  1  1 
1360  -0.1  1000  0.1  -1000 
1370  (ANCHORS) 

1380  C  6  90  2  1  1.E3 
1390  -0.043  800 
1400  C  17  70  2  1  1.E3 
1410  -0.052  10  0  0 
1420  FINISH 


Figure  36.  Input  file  for  Example  4 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  22-MAR-1994  TIME:  11.48.07 

*  INPUT  DATA  * 

I. — HEADING 

'EXAMPLE  4  —  MULTIPLE  ANCHORED  RETAINING  WALL 
'EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
'TENSION-ONLY  ANCHORS 


II.- 

-BEAM  CROSS 

SECTION  DATA 

< - 

- SECTION 

PROPERTIES 

X-COORDINATE 

MODULUS  OF 

< - 

START - > 

< - 

STOP - > 

START 

STOP 

ELASTICITY 

AREA 

INERTIA 

AREA 

INERTIA 

(IN) 

(IN) 

(PSD 

(SQIN) 

(IN**4) 

(SQIN) 

(IN**4) 

.  00 

552 . 00 

2.90E+07 

11.20 

280.80 

11.20 

280.80 

III. — NODE  SPACING  DATA 


X-COORDINATE 

MAXIMUM  NODE 

START 

STOP 

SPACING 

(FT) 

(FT) 

(FT) 

.  00 

50.00 

1.00 

IV. 

—  LOAD 

DATA 

NONE 

V.- 

-FIXED 

SUPPORT  DATA 

NONE 

VI. — LINEAR  SPRING  DATA 
NONE 


Figure  37.  Output  file  for  Example  4  (Sheet  1  of  8) 
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VII. — NONLINEAR  SPRING  DATA 


VII. A. — NONLINEAR  CONCENTRATED  SPRINGS 


SPRING 

X-COORD 

(FT) 

6.00 


ANGLE 

(DEG) 

90.00 


CURVE  COORD.  MULTIPLIERS 
DISPLACEMENT  FORCE 

(FT)  (P) 

l.OOE+00  l.OOE+03 


DISPLACEMENT 

COORDINATE 

-4.300E-02 

O.OOOE+00 


FORCE 

COORDINATE 

8.000E+00 

O.OOOE+00 


SPRING 

X-COORD 

(FT) 

17.00 


ANGLE 

(DEG) 

70.00 


CURVE  COORD.  MULTIPLIERS 
DISPLACEMENT  FORCE 

(FT)  (P) 

l.OOE+00  l.OOE+03 


DISPLACEMENT 

COORDINATE 

-5.200E-02 

O.OOOE+00 


FORCE 

COORDINATE 

l.OOOE+01 

O.OOOE+00 


VII. B. — NONLINEAR  DISTRIBUTED  SPRINGS 

SPRING  SPRING 

DISTRIBUTION  DIRECTION  X-COORD 

(FT) 

STARTS  Y  .00 


DISPLACEMENT 

COORDINATE 

-8.000E-03 

1.250E-01 


FORCE 

COORDINATE 

-1.330E+02 

-1.200E+03 


SPRING 

DISTRIBUTION  DIRECTION 


CONTINUES  Y 


SPRING 

X-COORD 

(FT) 

16.00 


DISPLACEMENT 

COORDINATE 

-8.000E-03 

1.250E-01 


FORCE 

COORDINATE 

-7.470E+02 

-6.720E+03 


CURVE  COORD. 
DISPLACEMENT 
(FT) 

l.OOE+00 


CURVE  COORD. 
DISPLACEMENT 
(FT) 

l.OOE+00 


Figure  37.  (Sheet  2  of  8) 


I 


MULTIPLIERS 

FORCE 

(P) 

l.OOE+00 


MULTIPLIERS 

FORCE 

(P) 

l.OOE+00 
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DISTRIBUTION 


CONTINUES 


SPRING 

DIRECTION 


DISPLACEMENT 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  22-MAR-1994  TIME:  11.48.24 


************************ 
*  SUMMARY  OF  RESULTS  * 
************************ 


I. — HEADING 

'EXAMPLE  4  —  MULTIPLE  ANCHORED  RETAINING  WALL 
'EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
'TENSION-ONLY  ANCHORS 


II.— MAXIMA 


AXIAL  DISPLACEMENT  (FT) 
LATERAL  DISPLACEMENT  (FT) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-FT) 


MAXIMUM 
POSITIVE 
8.584E-03 
1.775E-02 
6.585E-04 
3 . 128E+02 
7.412E+03 
3.618E+04 


X-COORD 

(FT) 

17.00 

.00 

46.00 

17.00 

17.00 

23.00 


MAXIMUM  X-COORD 

NEGATIVE  (FT) 

O.OOOE+00  .00 

-1.725E-01  32.00 

-8.925E-03  7.00 

-3.107E+03  17.00 

-4.155E+03  26.00 

-5.638E+03  6.00 


III. — FORCES  IN  CONCENTRATED  NONLINEAR  SPRINGS 


X-COORD 

ANGLE 

DEFORMATION 

(FT) 

(DEG) 

(FT) 

6.00 

90.00 

-3.452E-02 

17.00 

70.00 

-1.145E-01 

FORCE 

(P) 

-6.423E+03 

l.OOOE+04 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS 

WITH  NONLINEAR  SUPPORTS 

DATE: 

22-MAR-1994 

TIME:  11.48. 

24 

*  COMPLETE  RESULTS 

* 

I.— 

HEADING 

'EXAMPLE  4  —  MULTIPLE  ANCHORED  RETAINING  WALL 

'EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 

'TENSION-ONLY  ANCHORS 

II.— 

DISPLACEMENTS  AND  INTERNAL  FORCES 

X-COORD 

AXIAL  LATERAL  ROTATION 

AXIAL 

SHEAR 

MOMENT 

(FT) 

(FT)  (FT)  (RAD) 

(P) 

(P) 

(P-FT) 

.00 

8.578E-03  1.775E-02  -8.660E-03 

0.  OOOE+00 

O.OOOE+00 

O.OOOE+00 

1.00 

8.578E-03  9.087E-03  -8.661E-03 

5.986E+00 

-3.470E+02 

-1.728E+02 

2.00 

8.578E-03  4.236E-04  -8.668E-03 

1. 348E+01 

-6.824E+02 

-6.902E+02 

3.00 

8.578E-03  -8.253E-03  -8.687E-03 

2.248E+01 

-9.463E+02 

-1.510E+03 

4 . 00 

8.578E-03  -1.696E-02  -8.723E-03 

3.298E+01 

-1.214E+03 

-2.587E+03 

5.00 

8.578E-03  -2.571E-02  -8.780E-03 

4.499E+01 

-1.519E+03 

-3.950E+03 

6.00 

8.578E-03  -3.452E-02  -8.864E-03 

5.851E+01 

-1.863E+03 

-5.638E+03 

6.00 

8.578E-03  -3.452E-02  -8.864E-03 

5.851E+01 

4.560E+03 

-5.638E+03 

7.00 

8.579E-03  -4.343E-02  -8.925E-03 

7 . 354E+01 

4.177E+03 

-1.266E+03 

8.00 

8.579E-03  -5.235E-02  -8.912E-03 

9.007E+01 

3 .757E+03 

2.704E+03 

9.00 

8.579E-03  -6.123E-02  -8.832E-03 

1.081E+02 

3.297E+03 

6.234E+03 

10.00 

8.579E-03  -6.999E-02  -8.694E-03 

1.277E+02 

2.800E+03 

9.286E+03 

11.00 

8.580E-03  -7.860E-02  -8.506E-03 

1.487E+02 

2.264E+03 

1.182E+04 

12.00 

8.580E-03  -8.699E-02  -8.279E-03 

1.713E+02 

1.690E+03 

1.380E+04 

13.00 

8.581E-03  -9.515E-02  -8.022E-03 

1.953E+02 

1.077E+03 

1.519E+04 

14.00 

8.582E-03  -1.030E-01  -7.746E-03 

2.209E+02 

4 .259E+02 

1.594E+04 

15.00 

8.582E-03  -1.106E-01  -7.462E-03 

2.480E+02 

-2.635E+02 

1.603E+04 

16.00 

8.583E-03  -1.180E-01  -7.183E-03 

2.766E+02 

-9.913E+02 

1. 540E+04 

17.00 

8.584E-03  -1.250E-01  -6.922E-03 

3.128E+02 

-1.985E+03 

1. 395E+04 

17.00 

8.584E-03  -1.250E-01  -6.922E-03 

-3.107E+03 

7.412E+03 

1.395E+04 

18.00 

8.575E-03  -1.318E-01  -6.613E-03 

-3.064E+03 

6.189E+03 

2.075E+04 

19.00 

8.565E-03  -1.382E-01  -6.195E-03 

-3.019E+03 

4.962E+03 

2.633E+04 

20.00 

8.556E-03  -1.441E-01  -5.689E-03 

-2.972E+03 

3.725E+03 

3.067E+04 

21.00 

8.547E-03  -1.496E-01  -5.118E-03 

-2.924E+03 

2.475E+03 

3.378E+04 

22.00 

8.538E-03  -1.544E-01  -4.502E-03 

-2.875E+03 

1.207E+03 

3.562E+04 

23.00 

8.529E-03  -1.586E-01  -3.865E-03 

-2.824E+03 

-8.444E+01 

3.618E+04 

24.00 

8.521E-03  -1.621E-01  -3.230E-03 

-2.772E+03 

-1.405E+03 

3.544E+04 

25.00 

8.512E-03  -1.650E-01  -2.620E-03 

-2.718E+03 

-2.760E+03 

3.336E+04 

26.00 

8.504E-03  -1.674E-01  -2.058E-03 

-2.663E+03 

-4 . 155E+03 

2.991E+04 

27 . 00 

8.496E-03  -1.692E-01  -1.565E-03 

-2 . 563E+03 

-3.841E+03 

2.591E+04 

28.00 

8.488E-03  -1.705E-01  -1.140E-03 

-2.460E+03 

-3.520E+03 

2.223E+04 

29 . 00 

8.481E-03  -1.715E-01  -7.773E-04 

-2 . 353E+03 

-3.197E+03 

1.887E+04 

30.00 

8.474E-03  -1.721E-01  -4.710E-04 

-2 . 243E+03 

-2.876E+03 

1.583E+04 
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31.00 

8.467E-03  -1.724E-01 

-2 . 155E-04 

-2.129E+03 

-2 . 561E+03 

1. 311E+04 

32 . 00 

8.460E-03  -1.725E-01 

-5.342E-06 

-2.011E+03 

-2.255E+03 

1. 071E+04 

33.00 

8.454E-03  -1.724E-01 

1.649E-04 

-1.890E+03 

-1.962E+03 

8.599E+03 

34 . 00 

8.449E-03  -1.722E-01 

3 . 005E-04 

-1.766E+03 

-1. 684E+03 

6.777E+03 

35.00 

8.444E-03  -1.719E-01 

4 . 062E-04 

-1.638E+03 

-1.423E+03 

5.225E+03 

36.00 

8.439E-03  -1.714E-01 

4 . 867E-04 

-1.507E+03 

-1. 181E+03 

3.924E+03 

37 . 00 

8.434E-03  -1.709E-01 

5.464E-04 

-1.372E+03 

-9. 593E+02 

2.855E+03 

38 . 00 

8.430E-03  -1.703E-01 

5.890E-04 

-1.233E+03 

-7 . 587E+02 

1.998E+03 

39 . 00 

8.427E-03  -1.697E-01 

6. 182E-04 

-1.091E+03 

-5.803E+02 

1.331E+03 

40.00 

8.424E-03  -1.691E-01 

6.370E-04 

-9.461E+02 

-4 . 248E+02 

8 . 301E+02 

41.00 

8.421E-03  -1.684E-01 

6.484E-04 

-7.972E+02 

-2 . 928E+02 

4.733E+02 

42 . 00 

8.419E-03  -1.678E-01 

6. 545E-04 

-6.448E+02 

-1. 849E+02 

2 . 364E+02 

43.00 

8.417E-03  -1.671E-01 

6. 573E-04 

-4.889E+02 

-1. 013E+02 

9.538E+01 

44.00 

8.416E-03  -1.665E-01 

6. 583E-04 

-3.295E+02 

-4 . 244E+01 

2.557E+01 

45.00 

8.415E-03  -1.658E-01 

6.585E-04 

-1.665E+02 

-8.572E+00 

2 . 166E+00 

46.00 

8.415E-03  -1.652E-01 

6.585E-04 

O.OOOE+00 

0.  OOOE+00 

0. OOOE+00 

III. — FORCES  IN  DISTRIBUTED  NONLINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 


•COORD 

AXIAL 

LATERAL 

(FT) 

(P/FT) 

(P/FT) 

.00 

-5.233E+00 

-3.396E+02 

1.00 

-6.739E+00 

-3 .479E+02 

2.00 

-8.246E+00 

-3 . 162E+02 

3.00 

-9.752E+00 

-2.481E+02 

4 . 00 

-1.126E+01 

-2.865E■^02 

5.00 

-1.277E+01 

-3.249E+02 

6.00 

-1.427E-)-01 

-3.633E+02 

7.00 

-1.578E-f01 

-4 . 016E-(-02 

8 . 00 

-1.729E+01 

-4 .400E+02 

9.00 

-1.879E-f01 

-4 . 784E-t-02 

10.00 

-2.030E-t-01 

-5. 168E-f02 

11.00 

-2.181E+01 

-5.551E-t-02 

12.00 

-2.332E-t-01 

-5.935E-(-02 

13.00 

-2.483E+01 

-6. 319E-f02 

14.00 

-2.633E-f01 

-6.703E-f02 

15.00 

-2.784E+01 

-7.086E-f02 

16.00 

-2.935E-t-01 

-7.470E-t-02 

17.00 

-4.292E+01 

-1.221E+03 

18.00 

-4.435E+01 

-1.224E-t-03 

19.00 

-4.578E+01 

-1.231E+03 

20.00 

-4.721E-f01 

-1.242E-t-03 

21.00 

-4 . 863E+01 

-1.258E+03 

22.00 

-5.005E+01 

-1.279E-(-03 

23.00 

-5. 147E-f01 

-1.305E+03 

24 . 00 

-5.289E-f01 

-1.337E-t-03 

25.00 

-5.430E+01 

-1.374E+03 
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26.00 

-5.571E+01 

-1.416E+03 

26.00 

-9.823E+01 

3 . 081E+02 

27.00 

-1.017E+02 

3. 182E+02 

28.00 

-1.052E+02 

3.230E+02 

29.00 

-1.087E+02 

3.229E+02 

30.00 

-1. 122E+02 

3. 187E+02 

31.00 

-1.157E+02 

3. 108E+02 

32 . 00 

-1. 192E+02 

2 .997E+02 

33 . 00 

-1.227E+02 

2.859E+02 

34.00 

-1.261E+02 

2 . 698E+02 

35.00 

-1. 296E+02 

2.518E+02 

36.00 

-1. 331E+02 

2.323E+02 

37.00 

-1.366E+02 

2.115E+02 

38.00 

-1.401E+02 

1.897E+02 

39.00 

-1.436E+02 

1.671E+02 

40.00 

-1.471E+02 

1.438E+02 

41.00 

-1.506E+02 

1.200E+02 

42.00 

-1.541E+02 

9.582E+01 

43.00 

-1.577E+02 

7.128E+01 

44.00 

-1.612E+02 

4.642E+01 

45.00 

-1.647E+02 

2. 127E+01 

46.00 

-1.683E+02 

-4.171E+00 
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Figure  38.  Axial  and  bending  stiffnesses  for  Example  4 


•EXAMPLE  4  —  MULTIPLE  ANCHORED  RETAINING  WALL 
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NONLINEAR  CONCENTRATED  SPRING 


■EXAMPLE  4  --  MULTIPLE  ANCHORED  RETAINING  HALL 
■EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
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Rgure  40.  Schematic  of  nonlinear  distributed  springs  for  Example  4 


-  MULTIPLE  ANCHORED  RETAINING  HALL 
AXIAL  FORCE  ON  BENDING  EXCLUDED 
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Figure  41 .  Axial  and  lateral  displacements  for  Example  4 


■EXAMPLE  4  --  MULTIPLE  ANCHORED  RETAINING  HALL 
'EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 


Chapter  5  Example  Soultions 


Figure  42.  Shear  force  and  bending  moment  diagrams  for  Example  4 


■EXAMPLE  4  --  MULTIPLE  ANCHORED  RETAINING  WALL 
•EFFECTS  OF  AXIAL  FORCE  OH  BENDING  EXCLUDED 
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FORCES  IN  Y-DISTRIBUTED  NONLINEAR  SPRINGS  (UNITS:  P/FT) 
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Appendix  A 
Guide  for  Data  input 


Source  of  Input 

Input  data  may  be  supplied  from  a  predefined  data  file  or  from  the  user 
terminal  during  execution.  If  data  are  supplied  from  the  user  terminal,  prompt¬ 
ing  messages  are  printed  to  indicate  the  amount  and  character  of  data  to  be 
entered. 


Data  Editing 

When  all  data  for  a  problem  have  been  entered,  the  user  is  offered  the 
opportunity  to  review  an  echoprint  of  the  currently  available  input  data  and  to 
revise  any  or  all  sections  of  the  input  data  before  execution  is  attempted.  When 
editing  during  execution,  each  section  must  be  entered  in  its  entirety. 


Input  Data  File  Generation 

After  data  have  been  entered  from  the  terminal,  either  initially  or  after 
editing,  the  user  may  direct  the  program  to  write  the  input  to  a  permanent  file 
in  input  data  file  format. 


Data  Format 

All  input  data  (whether  supplied  fiom  the  user  terminal  or  from  a  file)  are 
read  in  free-field  format: 

a.  Data  items  must  be  separated  by  one  or  more  blanks  (comma  separators 
are  not  permitted). 

b.  Integer  numbers  must  be  of  the  form  NNNN. 
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c.  Real  numbers  may  be  of  the  form  ±xxxx,  ±xx.xx,  or  ±xx.xxE±ee. 


d.  User  responses  to  all  requests  for  control  by  the  program  for 

alphanumeric  input  may  be  abbreviated  by  the  first  letter  of  the  indi¬ 
cated  word  response,  e.g., 

ENTER  ’YES’  OR  ’NO’-respond  Y  or  N 

ENTER  ’CONTINUE’  OR  ’END’-respond  C  or  E. 


Sections  of  Input 

Input  data  are  divided  onto  the  following  sections; 

a.  Heading. 

b.  Beam  cross  section  properties. 

c.  Node  spacing. 

d.  Applied  loads. 

e.  Fixed  supports. 

/.  Linear  spring  supports. 

g.  Nonlinear  spring  supports. 

h.  Termination. 

When  data  are  entered  from  the  terminal,  data  sections  a  through  g  may  be 
input  in  any  order.  When  data  are  entered  from  a  predefined  data  file,  the 
heading  section  a  must  be  entered  first;  other  sections  b  through  g  may  be 
entered  in  any  order. 

When  data  are  entered  from  the  terminal,  the  user  is  cued  for  Termination. 
Termination  for  a  data  file  is  discussed  later. 


Minimum  Required  Data 

Data  sections  a,  b,  and  c  are  always  required.  At  least  one  of  sections  e,  /, 
and  g  is  required.  It  is  the  responsibility  of  the  user  to  ensure  that  sufficient 
supports  (fixed  or  spring  supports)  are  provided  to  inhibit  all  rigid  body 
motions  of  the  system,  i.e.,  to  ensure  a  stable  structure. 
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Units 


The  program  recognizes  the  following  units: 

/nches.  Feet,  Founds,  ATips 

Default  is  to  inches  and  pounds.  Angular  units  are  degrees  or  radians  as 
explained  below. 

Each  data  section  may  be  entered  with  any  combination  of  units  for  length 
and  force  as  specified  by  the  user. 

Predefined  Data  Fiie 

In  addition  to  the  general  format  requirements  given  in  paragraphs  entitled 
“Data  Format,”  the  following  pertain  to  a  predefined  data  file  and  to  the  input 
data  description  which  follows: 

a.  Each  line  must  commence  with  a  nonzero,  positive  line  number, 
denoted  LN  below. 

b.  A  line  of  input  may  require  both  alphanumeric  and  numeric  data  items. 
Alphanumeric  data  items  are  enclosed  in  single  quotes  in  the  following 
paragraphs. 

c.  A  line  of  input  may  require  a  keyword.  The  acceptable  abbreviation 
for  the  keyword  is  indicated  by  the  initial  italicized  capital  letters  of  the 
word.  The  acceptable  abbreviation  for  the  keyword  'UNeaf  is  'UN\ 
e.g. 

d.  Lower  case  words  in  single  quotes  indicate  a  choice  of  keywords 
defined  follows. 

e.  Items  designated  by  upper  case  letters  and  numbers  without  quotes 
indicate  numeric  data  values.  Numeric  data  values  are  either  real  or 
integer  according  to  standard  FORTRAN  variable  naming  conventions. 

/.  Data  items  enclosed  in  brackets  [  ]  may  not  be  required.  Data  items 
enclosed  in  braces  {  }  indicate  special  note  follows. 

g.  Input  data  are  divided  into  the  sections  discussed  in  paragraphs  entitled 
“Sections  of  Input”  Except  for  the  heading,  each  section  consists  of  a 
header  line  and  one  or  more  data  lines.  The  header  line  serves  the 
multiple  putposes  of:  indicating  the  end  of  the  preceding  section;  iden¬ 
tifying  the  data  section  to  follow;  indicating  whether  the  section  is 
’New’  or  to  be  ’Added’  to  the  same  section  of  a  preceding  problem 
(following  page);  and  indicating  the  units  associated  with  the  section. 
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h.  A  data  file  may  contain  input  data  for  several  problems  to  be  run  in 
succession.  On  the  first  problem  defined  in  the  file,  all  data  sections 
must  be  'New'.  In  subsequent  problem  descriptions,  'New'  indicates 
the  section  is  to  replace  that  section  in  the  preceding  problem;  ’Add’ 
indicates  the  data  are  to  be  appended  to  that  section  in  the  preceding 
problem.  A  data  section  may  be  deleted  from  the  preceding  problem 
description  by  indicating  'New'  on  the  section  header  line  and  immedi¬ 
ately  following  with  another  section  header. 

i.  Each  section  header  line  contains  the  information  [{’A^ew’  or  ’Add’}] 
and  ({’units’}].  Either  or  both  of  these  items  may  be  omitted. 

(1)  If  [{’A^ew’  or  ’Add’}]  is  omitted,  'New'  is  assumed.  Omit  this  item 
on  the  first  problem  description  in  the  file. 

(2)  If  ’Add’  is  indicated,  the  units  for  appended  data  must  be  the 
same  as  for  that  section  in  the  preceding  problem. 

(3)  [{’units’}]  must  be  'Inches'  or  ’Feet’  and  ’Founds’  or  ’ATips’.  If 
data  are  to  be  appended  to  the  preceding  problem,  omit  this  item. 
Data  to  be  appended  must  be  in  the  same  units  prescribed  for  that 
section  in  the  preceding  problem. 

(4)  If  {{’units’}]  is  omitted  in  the  header  for  a  'New'  section,  default 
is  to  inches  and  pounds. 

J.  If  data  are  added  to  a  preceding  section,  the  number  of  lines  in  the  pre¬ 
ceding  section  plus  the  number  of  lines  added  must  not  exceed  the 
limitations  on  total  lines  in  a  section  specified  in  the  input  description 
which  follows. 

k.  A  data  file  must  begin  with  a  Heading.  In  a  succession  of  problems, 
the  Heading  for  the  preceding  problem  may  be  replaced  by  entering  the 
new  heading  at  the  beginning  of  the  new  problem  description.  Addi¬ 
tion  to  a  previous  heading  is  not  pennitted. 

/.  Comment  lines  may  be  inserted  in  the  input  file  by  enclosing  the  line 
following  the  line  number  in  parentheses.  Comment  lines  are  ignored, 
e.g., 

1234  (THIS  LINE  IS  IGNORED) 


Input  Description 

Heading-One  (1)  to  four  (4)  lines  for  identifying  the  problem. 
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a.  Line  contents. 


LN  {’heading’} 

b.  Definition. 

’heading’  =  any  alphanumeric  infonnation  up  to  70  characters  including 
LN  and  any  imbedded  blanks. 

c.  Restriction:  If  a  ’heading’  line  following  LN  begins  with  a  combina¬ 
tion  of  letters  which  are  permissible  abbreviations  of  any  of  the  key¬ 
words  described  below,  Ae  ’heading’  must  begin  with  a  single  quote. 

Beam  Cross  Section  Properties-Two  (2)  to  twenty-two  (22)  lines. 

a.  Header-One  (1)  line. 

(1)  Contents. 

LN  ’Beam’  [{’New’  or  ’Add’}]  [{’units’}] 

(2)  Definition. 

’Beam’  =  section  title. 

[{’units’}]  =  ’/nches’  or  ’Feet’  and  ’Founds’  or  ’Nips’  for  ’New’; 
default  to  inches  and  pounds  if  omitted  for  ’New’. 

b.  Data  Lines— One  (1)  to  twenty-one  (21)  lines. 

(1)  Contents. 

LN  XI  X2  E  A1  SIl  [  A2  SI2] 

(2)  Definitions. 

XI  =  x-coordinate  at  beginning  of  distribution. 

X2  =  x-coordinate  at  end  of  distribution. 

E  =  modulus  of  elasticity  in  region  XI  to  X2. 

A1  =  cross  section  area  at  XL 
SIl  =  cross  section  moment  of  inertia  at  XL 
[A2,SI2]  =  area  and  moment  of  inertia  at  X2. 
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(3)  Discussion. 


(a)  If  A2  and  SI2  are  omitted,  area  and  moment  of  inertia  are 
assumed  to  be  constant  from  XI  to  X2. 

(b)  If  A2  and  SI2  are  provided,  area  and  moment  of  inertia  are 
assumed  to  vary  linearly  from  A1  and  SIl  at  XI  to  A2  and 
SI2  at  X2. 

(c)  When  several  distributions  are  described,  the  minimum  of  all 
XI  values  defines  the  left  end  of  the  beam  and  the  maximum 
of  all  X2  values  defines  the  right  end  of  the  beam.  There 
must  be  positive,  nonzero  values  of  modulus  of  elasticity, 
area,  and  moment  of  inertia  for  every  x-coordinate  between 
the  ends  of  the  beam. 

(d)  If  distributions  are  specified  for  overlapping  regions,  values 
within  the  overlap  are  cumulative. 


Node  Spacing-One  (1)  to  twenty-two  (22)  lines. 

a.  Header-One  (1)  line. 

(1)  Contents. 

LN  'NODc'  [{Wew’ or ’Add’}]  [{’units’}) 

(2)  Definitions. 

'NODt'  =  section  title. 

[{’units’}]  =  ’/nches’  or  ’Feet’  for  ’New’;  default  to  inches  if 
omitted  for  ’New’. 

b.  Data  Line-One  (1)  to  twenty-one  (21)  lines. 

(1)  Contents. 

LN  XI  X2  HMAX 

(2)  Definitions. 

XI, X2  =  x-coordinates  at  beginning  and  end  of  interval, 
respectively. 

HMAX  =  maximum  distance  between  nodes  in  interval  XI  to  X2. 


Appendix  A  Guide  for  Data  Input 


c.  Discussion. 


(1)  Intervals  XI  to  X2  must  proceed  sequentially  from  left  to  right. 
Overlapping  intervals  or  gaps  between  intervals  are  not  permitted. 

(2)  XI  on  the  first  line  input  must  be  less  than  or  equal  to  the 
x-coordinate  at  the  left  end  of  the  beam. 

(3)  X2  on  the  last  line  input  must  be  greater  than  or  equal  to  the 
x-coordinate  at  the  right  end  of  the  beam. 

(4)  Data  specified  beyond  the  ends  of  the  beam  are  ignored. 

(5)  Node  spacing  data  must  result  in  at  least  three  (3)  nodes,  (two 
(2)  elements)  and  not  more  than  two-hundred-one  (201)  nodes 
(200  elements)  in  the  model. 

Applied  Loads-Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines.  Entire  section 
may  be  omitted. 

a.  Header-One  (1)  line. 

(1)  Contents. 

LN  XOads’  [{Wew’ or  Add’}]  [{’units’}] 

(2)  Definitions. 

’LOads’  =  section  title. 

[{’units’}}  =  ’/nches’  or  ’Feet’  and  ’Founds’  or  ’Xips’  for  ’New’; 
default  to  inches  and  pounds  if  omitted. 

b.  Data  Lines  for  Concentrated  Loads-One  (1)  line  for  each  concentrated 
load. 

(1)  Contents. 

LN  ’Concentrated’  XI  FX  FY  C 

(2)  Definitions  for  concentrated  loads. 

’Concentrated’  =  keyword. 

XI  =  x-coordinate  at  point  of  concentrated  load(s). 

FXJ^  =  magnitudes  of  concentrated  x-  and  y-loads, 
respectively. 
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C  =  magnitude  of  concentrated  couple. 

c.  Data  Lines  for  Distributed  Loads-One  (1)  line  for  each  distribution. 

(1)  Contents. 

LN  ’Distributed’  {’direction’}  XI  Q1  X2  |Q2] 

(2)  Definitions. 

’Distributed  =  keyword. 

’direction’  =  ’X’  for  axial  load,  or  ’Y’  for  lateral  load. 

XI  =  x-coordinate  at  start  of  distribution. 

X2  =  x-coordinate  at  end  of  distribution. 

Q1  =  magnitude  of  distributed  load  at  XI. 

[Q2]  =  magnitude  of  distributed  load  at  X2:  assumed  to  be 
equal  to  Q1  if  omitted  (i.e.,  uniform  load). 

d.  Discussion. 

(1)  Multiple  loads  specified  at  a  single  x-coordinate,  either  concen¬ 
trated  loads  or  overlapping  distributed  loads,  are  cumulative. 

(2)  Data  specified  beyond  the  ends  of  the  beam  are  ignored. 

Fixed  Supports-Zeno  (0)  or  two  (2)  to  eleven  (11)  lines.  Entire  section 
may  be  omitted. 

a.  Header-One  (1)  line. 

(1)  Contents. 

LN  ’F/Xed’  [{’/Vew’ or ’/Idd’}]  [{’units’}] 

(2)  Definitions. 

'FlXtd'  =  section  title. 

{{’units’}]  =  ’/nches’  or  ’Feet’  for  ’A^ew’;  default  to  inches  if 
omitted. 

b.  Data  Lines-One  (1)  to  ten  (10)  lines. 

(1)  Contents. 
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LN  XI  {XD}  {YD}  {R} 


(2)  Definitions. 

XI  =  x-coordinate  at  point  of  support. 

(XD)  =  specified  X-displacement;  enter  Tree’  if  unspecified. 
(YD)  =  specified  Y-displacement;  enter  ’free’  if  unspecified. 

(R)  =  specified  rotation  (in  radians);  enter  ’free’  if  unspecified 
c.  Discussion. 

(1)  Multiple  values  of  a  single  displacement  component  are  not 
permitted. 

(2)  Data  specified  beyond  the  ends  of  the  beam  are  ignored. 

Linear  Spring  Supports--Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines. 
Entire  section  may  be  omitted. 

a.  Header-One  (1)  line. 

(1)  Contents. 

LN  ’L/Near’  [{’New’ or ’>ldd’}]  [{’units’}] 

(2)  Definitions. 

'UNeaf  =  section  title. 

[{’units’}]  =  ’/nches’  or  ’feet’  and  ’founds’  or  ’Xips’  for  ’New’; 
default  to  inches  and  pounds  if  omitted. 

b.  Data  Lines  for  Concentrated  Linear  Springs—One  (1)  line  for  each 
concentrated  linear  spring. 

(1)  Contents. 

LN  ’Concentrated’  XI  ANGLE  ST  SR 

(2)  Definitions. 

’Concentrated’  =  keyword. 

XI  =  x-coordinate  at  point  of  attachment. 
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ANGLE  =  angle  (in  degrees)  between  x-axis  and  line  of  action  of 
translation  spring;  measured  positive  counter-clockwise 
from  x-axis. 

ST  =  stiffness  of  translation  spring  in  force  per  unit 
displacement. 

SR  =  stiffness  of  rotation  spring  in  moment  per  radian. 

c.  Data  Lines  for  Distributed  Linear  Springs— One  (1)  line  for  each 
distribution. 

(1)  Contents. 

LN  ’Distributed’  {’direction’}  XI  X2  A  B  C 

(2)  Definitions. 

’Distributed’  =  keyword. 

{’direction’}  =  ’X’  for  axial  spring  or  ’Y’  for  lateral  spring. 

XI  =  x-coordinate  at  start  of  distribution. 

X2  =  x-coordinate  at  end  of  distribution. 

A.B.C  =  spring  stiffness  coefficients  (see  following 
discussion). 


d.  Discussion. 

(1)  If  several  concentrated  linear  springs  are  attached  to  a  single  point, 
effects  are  cumulative. 

(2)  Stiffness  of  a  distributed  spring  is  assumed  to  vary  according  to 

S  =  A  -H  B  •  tP,  where  Z  is  the  distance  from  XL  Units  of  A  are 
force  per  unit  length  per  unit  displacement;  units  of  B  and  C  must 
result  in  the  units  of  force  per  unit  length  per  unit  displacement  of 
the  product  B  •  Z^. 

(3)  Overlapping  distributions  result  in  cumulative  effects  in  the 
overlap. 

(4)  Data  specified  beyond  the  end  of  the  beam  are  ignored. 

Nonlinear  Spring  Supports-Zero  (0)  lines  or  maximum  of  forty  one  (41) 
nonlinear  concentrated  and/or  distributed  support  curves;  entire  section  may  be 
omitted. 
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a.  Control-One  (1)  line. 

LN  WOMinear’  [{Wew’  or  Mdd’}]  [{’units’}] 

b.  Definitions. 

’A^OMinear’  =  section  title. 

[{’units’}]  =  ’/nches’  or  ’Feet’  and  ’Founds’  or  ’ATips’;  default  to 
inches  and  pounds  if  omitted. 

c.  Data  Lines  for  Concentrated  Nonlinear  Springs-Two  (2)  or  more  lines 

for  each  concentrated  spring. 

(1)  Line  1  Contents: 

LN  ’Concentrated’  XI  ANGLE  NPTS  DMUL  FMUL 
Line  2  Contents; 

LN  DEF(1)  FORCE(l)  DEF(2)  FORCE(2)  [.  .  .  DEF(NPTS) 
FORCE(NPTS)] 

(Pairs  of  DEF( )  FORCE( )  may  be  continued  on  succeeding  cards 
following  a  line  number  until  NPTS  pairs  have  been  entered.) 

(2)  Definitions. 

’Concentrated’  =  keyword. 

XI  =  x-coordinate  at  point  of  attachment. 

ANGLE  =  angle  (in  degrees)  between  x-axis  and  spring  line 
of  action  measured  positive  counterclockwise 
from  positive  x-direction. 

NPTS  =  number  of  points  on  resisting  force-deformation 
curve;  minimum  of  two  (2)  points  required;  max¬ 
imum  of  twenty  one  (21)  points  permitted. 

DMUL  =  scale  factor  for  curve  deformation  coordinate 
values;  must  be  positive  and  nonzero. 

FMUL  =  scale  factor  for  curve  resisting  force  coordinate 
values. 

DEF( )  =  curve  deformation  coordinate  value;  must  pro¬ 
ceed  sequentially  from  most  negative  to  most 
positive. 
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FORCE( )  =  curve  resisting  force  coordinate  value  corresponding 
to  DEF( ). 

d.  Data  Lines  for  Distribution  Nonlinear  Springs-Minimum  of  Two  (2) 

groups  for  each  distribution. 

(1)  Group  1,  Line  1  Contents: 

LN  ’Distributed’  {’direction’}  XI  NPTS  DMUL  FMUL 

Group  1,  Line  2  Contents: 

LN  DEF(l)  FORCE(l)  DEF(2)  FORCE(2)  [.  .  .  DEF(NPTS) 

FORCE(NPTS)] 

(Pairs  of  DEF( )  FORCE( )  may  be  continued  on  succeeding  cards 

following  a  line  number  until  NPTS  pairs  have  been  entered.) 

(2)  Definitions  for  Group  1: 

’Distributed’  =  keyword  indicating  beginning  of  distribution. 

{’direction’}  =  'X'  or  ’F  indicating  distributed  spring  resists  axial 
(x)  or  lateral  (y)  displacements. 

XI  =  x-coordinate  at  beginning  of  distribution. 

NPTS  =  number  of  points  provided  on  nonlinear  resisting 
force-displacement  curve  at  XI;  minimum  of  two 
(2)  points  required;  maximum  of  twenty  one  (21) 
points  permitted. 

DMUL  =  scale  factor  for  curve  deformation  coordinate  val¬ 
ues;  positive,  nonzero. 

FMUL  =  scale  factor  for  curve  resisting  force  coordinate 
values. 

DEF( )  =  deformation  coordinate  value;  must  proceed  from 
most  negative  to  most  positive. 

FORCE( )  =  resisting  distributed  force  corresponding  to  DEF( ). 

(3)  Group  2,  Line  1  Contents: 

LN  {’control’}  X2  NPTS  DMUL  FMUL 

Group  2,  Line  2  Contents: 
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LN  DEF(1)  FORCE(l)  DEF(2)  FORCE(2)  [.  .  .  DEF(NPTS) 
FORCE(NPTS)] 

(Pairs  of  DEF( )  Force( )  may  be  continued  on  succeeding  cards 
following  a  line  number  until  NPTS  pairs  have  been  entered.) 

(4)  Definitions  for  Group  2: 

{’control’}  =  ’Continue’  or  ’End’;  ’Continue’  indicates  distribu¬ 
tion  continues  and  another  set  of  Group  2  lines  fol¬ 
lows;  ’End’  indicates  this  is  the  last  curve  in  this 
distribution. 

X2  =  x-coordinate  for  this  nonlinear  curve. 

NPTS,  DMUL,  FMUL,  DEF(  ),  FORCE(  )  as  in  Group  1 
c.  Discussion. 

(1)  Final  deformation  and  resisting  force  curve  coordinates  are  prod¬ 
ucts  DMUL-DEF( )  and  FMlJL-FORCE( ),  respectively. 

(2)  Resisting  force-deformation  curve  must  be  single  valued  at  every 
point 

(3)  Resisting  force-deformation  curve  is  assumed  to  vary  linearly 
between  input  points. 

(4)  Resisting  force  is  assumed  to  be  constant  at  first  or  last  value  pro¬ 
vided  on  curve  for  deformations  beyond  the  range  of  deformation 
coordinates  provided. 

(5)  Several  concentrated  nonlinear  springs  may  be  attached  at  a  single 
point 

(6)  Characteristics  of  distributed  nonlinear  springs  at  intermediate 
points  on  a  distribution  are  obtained  by  linear  interpolation 
between  adjacent  curves. 

(7)  Overlapping  distributions  result  in  cumulative  effects  in  the 
overlap. 

(8)  Data  specified  beyond  the  ends  of  the  beam  are  ignored. 
Termination-One  (1)  line. 

a.  Contents. 

LN  ’E/Msh  [{’Eerun’}  [{’Eeep’}]] 
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b.  Definitions. 


T/ATish’  =  keyword  indicating  end  of  input  data  for  a  problem. 
Initiates  data  checking  and  solution. 

[{’7?erun’}]  =  keyword;  if  omitted,  program  assumes  this  is  end  of  data 
file;  if  included,  after  solution  of  preceding  problem, 
program  will  immediately  read  succeeding  data  as  input 
for  a  new  problem. 

[{’ATeep’}]  =  keyword;  if  omitted,  program  begins  with  zero  displace¬ 
ments  for  initial  iteration  of  nonlinear  problem;  if 
included,  iterations  begin  with  displacements  from  pre¬ 
ceding  solution  in  this  sequence  of  problems;  omit  if 
( { ’R^erun’ }  ]  is  omitted. 


Abbreviated  Input  Guide 

(Data  items  enclosed  in  brackets  [  ]  may  be  omitted.  Braces  {  }  enclosing 
data  lists  indicate  choose  one;  arrow  indicates  default  if  item  is  omitted.) 

Heading-One  (1)  to  four  (4)  lines. 

LN  ’heading’ 

[LN  ’heading’] 

[LN  ’heading’] 

[LN  ’heading’] 

Beam  Cross  Section  Properties-Two  (2)  to  twenty-two  (22)  lines. 
a.  Header-One  (1)  line. 


LN  'Beam' 


■^'New' 

'Add' 


-^'Inches' 

'het' 


I  -^'Pounds'} 

1  J 


b.  Data  Lines-One  (1)  to  twenty-one  (21)  lines. 

LN  XI  X2  E  A1  SIl  [A2  SI2] 

Node  Spacing— Two  (2)  to  twenty-two  (22)  lines. 
a.  Header-One  (1)  line. 
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.  r 1  r \^'Inches'\  ] 


b.  Data  Lines-One  (1)  to  twenty-one  (21)  lines. 

LN  XI  X2  HMAX 

Applied  Loads-Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines. 
a.  Header-One  (1)  line. 


LN  'leads' 


■^Wew’Y 
'Add'  I 


>'Inches']  j-^'^ounds' 
'Feet'  [  1  'Feet' 


b.  Data  Lines  for  Concentrated  Loads-One  (1)  line  for  each  concentrated 
I  load. 

I 


LN  ’Concentrated’  XI  FX  FY  C 


c.  Data  Lines  for  Distributed  Loads— One  (1)  line  for  each  distribution. 


LN  'Distributed' 


|’,^j  XI  Q1  X2  [Q2] 


Fixed  Supports-Zero  (0)  or  two  (2)  to  eleven  (11)  lines. 
a.  Header-One  (1)  line. 


b.  Data  Lines-One  (1)  to  ten  (10)  lines. 


LN  «  j  {tL] 


Linear  Spring  Supports-Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines. 
a.  Header-One  (1)  line. 


LN  'UNear' 


->’Nw’Y 
'Add'  J 


>'lnches'\  ^'^unds'\ 
'het'  j|  'Feet'  J 
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b.  Data  Lines  for  Concentrated  Linear  Springs-One  (1)  line  for  each 
spring. 

LN  ’Concentrated’  XI  ANGLE  ST  SR 

c.  Data  Lines  for  Distributed  Linear  Springs-One  (1)  line  for  each 
distribution. 


LN  'Distributed' 


XI  X2  A  B  C 


Nonlinear  Spring  Supports— Zero  (0)  or  maximum  of  forty  one  (41) 
nonlinear  concentrated  and/or  distributed  support  curves. 

a.  Header-One  (1)  line. 


b.  Data  Lines  for  Concentrated  Nonlinear  Springs-Two  (2)  or 
more  lines  for  each  spring. 

Line  1: 

LN  ’Concentrated’  XI  ANGLE  NPTS  DMUL  FMUL 
Line  2; 

LN  DEF(l)  FORCE(l)  DEF(2)  FORCE(2)  [.  .  .  DEF(NPTS) 
FORCE(NPTS)] 

c.  Data  Lines  for  Distributed  Nonlinear  Springs-Minimum  of  two 
(2)  groups  lines  for  each  distribution. 

Group  l.Line  1: 


LN  'distributed' 


XI  NPTS  DMUL  FMUL 


Group  1,  Line  2: 

LN  DEF(1)  FORCE(l)  DEF(2)  FORCE(2)  [.  .  .  DEF(NPTS) 
FORCE(NPTS)] 

Group  2,  Line  1 : 
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LN 


Continue']^ 

[  '§!^'  J 


X2  NPTS  DMUL  FMUL 


Group  2,  Line  2: 

LN  DEF(1)  FORCE(l)  DEF(2)  FORCE(2)  [.  .  .  DEF(NPTS) 
FORCE(NPTS)] 


Termination-One  (1)  line. 

LN  ’Finish’  [’Rerun’  [’Reep’}} 
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